Optical Sum Rule in Finite Bands 
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In a single finite electronic band the total optical spectral weight or optical sum carries information 
on the interactions involved between the charge carriers as well as on their band structure. It varies 
with temperature as well as with impurity scattering. The single band optical sum also bears some 
relationship to the charge carrier kinetic energy and, thus, can potentially provide useful information, 
particularly on its change as the charge carriers go from normal to superconducting state. Here we 
review the considerable advances that have recently been made in the context of high T c oxides, 
both theoretical and experimental. 
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I. INTRODUCTION 

Much is known about the properties of the supercon- 
ducting state in the cuprates, yet after 20 years of intense 
research there is, as yet, no consensus as to the driving in- 
teractions responsible for the pairing (mechanism). It is 
well established that the superconducting order param- 
eter has (i-wave rather than conventional s-wave sym- 
metry. Angular resolved photoemission spectroscopji^ 
(ARPES) finds that the leading edge of the electron spec- 
tral density at the Fermi energy shifts below the chemical 
potential as a result of the onset of superconductivity. 
The observed shift is zero in the nodal direction (tt, tt) 
and maximum in the antinodal direction (0, tt). More re- 
cently a non monotonic increase of the gap amplitude as 
a function of direction has been observed in an electron 
doped cuprate but the symmetry is again <i-wave£ 

Many other independent experimental techniques have 
provided additional compelling evidence for d-wave sym- 
metry. One such technique is microwave spectroscopy 
which gives a linear in temperature (T) dependence of 
the superfluid density at low T £| . Complimentary to the 
above techniques are phase sensitive experiments- based 
on Josephson tunneling and flux quantization which pro- 
vide direct evidence for a change in sign of the order 
parameter not just that it has a zero. 

Another important observation which points directly 
to a non conventional mechanism is the so called collapse 
of the inelastic scattering rateSiZ&S in the superconduct- 
ing state which results in a large peak at intermediate T 
below the critical temperature (T c ) in the real part of the 
microwave conductivityALifiiii The origin of this peak is 
understood in terms of a readjustment in the excitation 
spectrum involved in the interaction as superconductivity 
develops. This is expected in electronic mechanisms 6,9 
involving the spin or charge susceptibility which looses 
weight at small energies corresponding to a hardening 
of the spectrum and little inelastic scattering remains at 
low T. 

All the evidence reviewed above points to deviations 



from a simple BCS s-wave phonon driven superconduc- 
tivity, but, nevertheless, the search goes on for other es- 
sential differences which could help identify the under- 
lying mechanism. An avenue to explore is the idea of 
kinetic as opposed to potential energy driven supercon- 
ductivity. The idea is that it is a reduction in kinetic 
energy (KE) that drives the condensation into the su- 
perconducting state in contrast to BCS theory where the 
potential is decreased sufficiently to overcome a KE in- 
crease and provides as well a reduction in total energy. 
The possibility of KE driven superconductivity was con- 
sidered by Hirschi^ and explicitely demonstrated for the 
hole mechanism of superconductivitj*iii^*i^ for which the 
effective mass of the holes decreases in the superconduct- 
ing state by pairing. It is also expected in other theoret- 
ical frameworks for highly correlated systems>i£ii£ 

In a single finite tight binding band with nearest neigh- 
bors hopping only, the KE is related to the optical sum 
defined as the total optical spectral weight under the real 
part of the optical conductivity ai(T,u>) integrated over 
energy uji& This fact applies whatever the interactions 
involved, the temperature, and also when impurities are 
present. Of course, for this to hold the conductivity needs 
to be computed exactly including vertex corrections. At 
first sight the restriction to a tight binding band with 
nearest neighbors only appears restrictive but in reality 
it has been found that KE and optical sum (OS) track 
each other reasonably well even when higher neighbors 
are considered (second, third, etc.). In view of this fact 
several recent experimental and theoretical papers have 
appeared concerning the temperature dependence of the 
OS in the normal state and its change in the supercon- 
ducting state. Here, we review this body of work with 
the aim of understanding what it tells us about the un- 
derlying interactions involved. We will restrict our dis- 
cussion to the cuprates and to their in-plane response. 
For the c-axis motion a sum rule violation was noticed 
by Basov et alm^ but the KE involved in the c-axis mo- 
tion is smalli22*2i To investigate its relation to the total 
condensation energy of the transition to the supercon- 
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ducting state, it is necessary to consider the a6-plane 
of the cuprates. Exactly how to treat the c-axis trans- 
port also complicates the interpretation. Some possi- 
ble mechanisms include strong intra layer scattering, 22 
non Fermi liquid ground states^ confinement^ inter- 
plane and in-plane charge fluctuations fSa indirect o 
axis coupling through the particle-particle channel, 25 
resonant tunneling on localized states in the blocking 
layer j— two band models ^L2& and coherent and incoher- 
ent tunneling. 2 Si22*2ii224^ Clearly this problem is worth 
further study but is not part of this review which will 
deal only with the cuprates and their afr-plane response. 

The in-plane optical sum (OS) in the cuprates is ob- 
served to increase as the temperature is decreased below 
room temperature. In the normal phase it is often, but 
not always, found to follow a T 2 law. The changes are 
small but larger than is expected on a rigid band model 
neglecting interactions. In the superconducting phase, on 
the other hand, a change in slope of this T 2 behavior is 
seen. While there remain some differences in details be- 
tween the various experimental groups, for underdoped 
samples it is agreed that the OS falls above the extrapo- 
lated normal state value. This is clearly seen in Fig. PHI 
which was reproduced from van der Marel et alm^ and 
deals with Bi 2 Sr 2 CaCu 2 8+a; (BSCCO, Bi2212). This 
behavior is opposite to that expected on the basis of BCS 
theory for which the OS is predicted to decrease. This 
fact can be traced to an increase in KE due to the for- 
mation of Cooper pairs. On the other hand, such a con- 
ventional behavior of the OS or change in KE on enter- 
ing the superconducting state was observed in overdoped 
samples of BSCCO by Dcutscher et alm^ Their data is 
reproduced in our Fig. [2U] These authors also find that 
the crossover from negative to positive KE change occurs 
around, but slightly above optimum doping (see Fig. l21f> . 
These findings have been largely confirmed by Carbone et 
alm& It is these facts that this review is directly concerned 
with and seeks to understand. Theoretically we will find 
that the temperature dependence of the OS can in some 
cases be dominated by a term which is proportional to a 
specific average over energy of the quasiparticle inelastic 
scattering rates and this need not give a T 2 law. This 
term is missing in all theories that do not explicitly treat 
damping effects. Further, this non T 2 dependence has 
implications for the accuracy of experimental results on 
the OS difference between superconducting and normal 
state which require an extrapolation to low temperatures 
of normal state data known only above T c . 

The paper is structured as follows. We begin with 
theoretical considerations and then review experimental 
information as it relates to the calculations. Section Hi Al 
introduces the optical sum and its relation to the kinetic 
energy. In Sec. Ill Bl we provide simple analytic formulas 
for the KE and OS in a free electron model which allows 
one to understand some qualitative aspects of this rela- 
tionship. It also contains a general formulation of the 
OS for a simplified tight binding model which averages 
over anisotropies and greatly simplifies the mathematics. 



We argue that the temperature dependence of the OS is 
governed mainly by the inelastic scattering. Both, real 
and imaginary part of the charge carrier self energy are 
important. In Sec. IHCl we describe the Nearly Antifcr- 
romagnetic Fermi Liquid model (NAFFL), give the set 
of generalized Eliashberg equations needed for numeri- 
cal work, and also specify the electron-spin fluctuations 
interaction (MMP model). We give results for the tem- 
perature dependence of the OS and of the KE in the 
normal and superconducting state. In Sec. Ill PI we in- 
troduce the Hubbard model and Dynamical Mean Field 
Theory (DMFT) and give results. 

In Sec. IIII Al we summarize some of the important ef- 
fects a finite band cutoff has on the self energy. Sec- 
tion ^^B] is devoted to results for the OS when a delta 
function, and Sec. IIII CI when an MMP form (spin fluc- 
tuations) is used for the electron-boson interaction. Sec- 
tion [^^D] provides an analysis of the temperature depen- 
dence of the OS for coupling to a low energy boson which 
results in a linear rather than quadratic in T law. It is 
argued that other temperature dependences could arise 
when different model interactions are used. 

In Sec. IIV Al we investigate within the NAFFL model 
the effect on the OS of the collapse of the inelastic scat- 
tering on entering the superconducting state which can 
lead to an increase in the OS rather than the usual BCS 
behavior. We also present experimental results in the 
BSCCO cuprates. Section IIVBI deals with a simpli- 
fied qualitative model based on a temperature dependent 
scattering rate which decreases with T as T 4 . While this 
model is not accurate, it shows clearly how the KE in 
the superconducting state can decrease below its normal 
state value due to the collapse of the inelastic scattering. 
In Sec. IIV O we describe a related, more phenomenolog- 
ical model due to Norman and Pepin which has several 
common elements with the model discussed in Sec. IIV Al 
We also provide comparison with experiment and addi- 
tional theoretical results for the KE change on conden- 
sation into the superconducting state in Sec. lIVDl Sec- 
tion [^H^ gives DMFT results for the normal as well as 
superconducting state in the t-J model. KE and poten- 
tial energy are discussed and compared with experiment. 
Further results based on the Hubbard model are com- 
mented on as are those based on the negative U Hubbard 
model used to describe the BCS - BE (Bose - Einstein) 
crossover. 

Section |V] deals with models of the pseudogap state 
above T c (the superconducting critical temperature) that 
exists up to a temperature T* > T c , with T* the pseudo- 
gap temperature. In Sec. IV Al we discuss KE changes in 
the preformed pair model in which pairs form at the pseu- 
dogap temperature T* and superconductivity sets in only 
at a lower temperature when phase coherence is estab- 
lished. In Sec. IVBl we consider an alternative model to 
phase fluctuations of the pseudogap state, namely the D- 
density wave model, a competing ordered state with bond 
currents and associated magnetic moments. In Sec. IVII 
we deal briefly with the problem of spectral weight dis- 
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TABLE I: The two tight binding models used within this pa- 
per. Model A corresponds to the tight binding model dis- 
cussed by van der Marel et alw— t and t' are given in meV, 
the critical temperature T c in K, and the filling (n) is defined 
in Eq. 



Model 


t 


t' 


(n) 


T c 


A 


148.8 


40.9 


0.425 


90 


B 


100.0 


16.0 


0.4 


100 



tribution as a function of energy. A short summary is 
provided in Sec. IVIII 



II. THEORY 



General considerations 



The single band OS is defined a: 
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Here <Ji xx (lo) is the real part of the (x, x)-component 
of the optical conductivity, e is the electron charge, Ti 
Planck's constant, k momentum, a spin, and V the crys- 
tal volume. The limits — fl to f2 on the conductivity 
integral are to be taken to include all the contributions 
to <Ji xx (uj) from the band of interest and from no other. 
The integration from — Vt to fi can be restricted to (0, fi) 
by making use of the fact that the real part of the con- 
ductivity is an even function of to. Often the integra- 
tion on f2 is also extended to infinity under the condition 
that only transitions from the one band of interest are 
included. The sum over k extends over the entire first 
Brillouin zone, n^ a is the probability of occupation of 
a state |k, ex), and £k is the electron dispersion relation. 
W is the OS divided by ne 2 /h 2 and will be quoted in 
meV. In deriving this OS rule the one body part of the 
Hamiltonian gives the single particle orbitals of energy ek 
and the two body piece provides correlation effects be- 
yond those included in £k (referred to as kinetic energy 
even though it can contain effects of interactions). If 
all bands are included rather than just the one in Eq. Q 
one would obtain the exact sum rule proportional to n/m 
with m the bare electron mass and n the total electron 
density. This is a constant independent of temperature. 
It does not change when the system undergoes a phase 
transition to the superconducting state by charge con- 
servation. This also applies to the free electron infinite 
band and is the basis for the Ferrell-Glover-Tinkham sum 
rul o 37 i 38 which will be elaborated upon in Sec. IVII 

The important observation for this review is that for 
tight binding bands with nearest neighbors hopping, i, 
only, i.e.: 
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FIG. 1: Optical sum W and kinetic energy, — VKkb/2, as a 
function of T 2 . Solid circles and squares are for the non inter- 
acting case while solid up-triangles and solid down-triangles 
include interactions. The left hand frame applies to Model A 
of Table U] and an MMP model parameter ljsf = 82meV [see 
Eq. 119H 1 was used. Here, the interacting and non interact- 
ing cases show similar temperature dependencies. The right 
hand frame is for Model B of Table Q with an MMP model 
parameter ljsf = lOmeV. Note the difference in temperature 
dependence between interacting and non interacting case. 



where a is the Cu02 plane lattice parameter, 
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1, 



-W KE 
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(3) 



Wke is the kinetic energy per copper atom. Thus, for 
this simple case an experimental measurement of W gives 
a measure of the KE. The band structure in the oxides 
is not, in general, limited to first neighbor hopping and 
one may well wonder if the above observation is of any 
practical use. 

For second neighbor hopping (i, t' model) 



£k = —2 {t [cos(ak x ) + cos(aky)] 
—It' cos(ak x ) cos(afcy)} — /i, 



(4) 



-2t [cos(afca;) + cos(afcj,)] 



(2) 



with n the chemical potential. In Fig. ^ we compare 
results for W (solid circles) with results for — Wke/2 
(solid squares) for two models, A and B, left and right 
hand frame, respectively. Parameters of the models arc 
defined in Table [I] where (n) is the filling factor which 
determines the chemical potential; with (n) = 0.5, half 
filling, /i = and t' = in Eq. J3J. No interactions are 
included beyond those giving rise to the tight binding 
band structure. It is clear that while W in both cases is 
a few percent smaller than — Wke/2 their dependence on 
temperature track each other fairly well. Thus, it makes 
sense to pursue this avenue as a means to get informa- 
tion on the KE and its variation with temperature and 
also when it undergoes a phase transition to the super- 
conducting state. One should be aware, however, that 
for some sets of tight binding parameters the difference 
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between OS and KE can be much larger and they may 
no longer track each other (see the very recent paper 
by Marsiglio et al£^ for more discussion). A final point 
should be stressed. To compute the OS in Eq. it is 
not necessary, although one can if one wishes, to calculate 
the conductivity which depends on a two-particle Green's 
function and is to be calculated with vertex corrections. 
The right hand side of Eq. Q , instead, depends only on 
the one particle spectral density which determines 
the probability of occupation of the state |k, a). On the 
other hand, if one is interested in the optical spectral 
weight integrated only to lom with ujm < one can no 
longer avoid calculating the conductivity. For this reason 
the problem of optical weight distribution is left to a last 
section IvTl 



B. The Continuum Limit and Quadratic 
Dispersion; Simple Results 

To get some insight into the significance of the OS, we 
consider next several simplifications. The probability of 
occupation of the state |k) (both spins) is related to the 
charge carrier spectral density A(k.,u>) by 



n k (T) 



J du f(T,u)A{k,uj), 



(5) 



where f(T, w) is the Fermi-Dirac distribution function 
at temperature T . For tight binding electrons with no 
correlation effects the spectral function ^4(k, u>) = <5(ek — 
w) and rik(T) = /(ek,T). Before proceeding to include 
interactions it is convenient in what follows to rewrite 
Eq. in an equivalent form through integration by parts 
with the surface integral equal to zero by symmetry in a 
periodic band: 



W = - 



k 



<9k 



drik 



(6) 



Knigavko et a/. 40 have found this second form to be 
very useful when making approximations. In the Kubo 
formula which defines the conductivity cri xx (u>) which ap- 
pears on the left hand side of Eq. (|TJ it is the velocity 
squared, (de^/dk x ) 2 that enters naturally. In making ap- 
proximations to the band structure in order to get sim- 
plified expressions it is important to have the same quan- 
tity (5gk/ 'dk x ) 2 on both sides of Eq. (JJJ and, indeed, in 
Ref. |4(J it is verified that the right hand side (RHS) of 
Eq. (JTJ and W of Eq. © agree to high accuracy for the 
two simplified band structure models they considered, 
namely free electron bands cut off to —D/2 and D/2 at 
half filling and another model designed to treat the tight 
binding case which we will define later. In both cases the 
density of states was approximated by a constant. Such 
an approximation, however, is not completely compatible 
with the assumed Fermi velocity model and this explains 



why we use W rather than W . In the continuum limit of 
Eq. (J2J with t' = fi = as the lattice parameter a — > 
one gets 

f-2 



£k 



-At 



2m 



{k x + k y ) , 



with Ti 2 /2m = ta 2 (where this product is assumed to be 
finite as a — > 0) or a 2 — — 4 and D/2 = At. For these 

> m D I 

approximations 

m v ) 2 



w = w 



Aire 2 



the well known result for free electron bands, with U 2 = 
Anne 2 jm the plasma frequency squared and n the elec- 
tron density. We can also show that 

W KE = -W, 

so that the relationship between KE and OS found for 
tight binding bands with first neighbors is profoundly 
changed when these are approximated by their contin- 
uum limit (free electron case). This is, in a sense, the 
extreme opposite case to the nearest neighbor only tight 
binding model. These results were obtained at zero tem- 
perature. At finite T things are not quite as simple. W 
and W undergo no change with temperature but Wke 
does and is 



W KE = 



N(0)Dh 2 



tt 2 fk B TY 
3 \D/2) 
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which is the well known result that the KE increases 
as the square of the temperature normalized to the half 
band width D/2. These results show that the correspon- 
dence between W and Wke can be rather subtle and 
it can be lost when approximations to the tight binding 
band are made. Here N(0) is the charge carrier density 
of states at the Fermi energy. 

Knigavko et aL— used a somewhat more sophisticated 
model to approximate a tight binding band. This model 
was used previously by Marsiglio and Hirsch^i In this 
approach the square of the electron velocity is replaced 
by 
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J 



Here, m b is a band electronic mass. In this model the rule 
that W = -Wke/2 still holds. For A(e,u>) = S(e - a) 
(no correlations) we recover 



W(T) = 



m b 



3 \D/2) 



(8) 



which is a result obtained by van der Marel et al& and 
often used to interpret experiments. It is important to 
contrast this result with the free electron case for which 
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duf(T,u)A(D/2,w) 



(9) 



where n is the charge carrier density. 

To treat interactions it is convenient to rewrite Eq. Q 
in the form 



W(T) 



2m(, 



where 



h(T,u) 



(D/2)' 



dw f(T,w)h(T,u) 



D/2 

J de eA(e,uj). 



(10) 



(11) 



Note also that the probability of occupation of the state 
e given in Eq. (J5J is closely related to Eqs. ljTU|) and (tTTt . 
We can apply the Sommerfeld expansion to Eq. Ijllfl to 
get 



W(T) 



2mj 

dh(T,u>) 



1 



dujh(T,u)-^-{k B T) 2 
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duo 
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(12) 



For no interactions h(T,u) — Alo/ (D /2) 2 for u> > and 
zero for u> < [note that the derivative in Eq. I|12l) (non 
interacting case) is to be weighted by 1/2 at u> = 0] so 
that the second term in Eq. I|12(l vanishes and the third 
term gives the first correction in |JS}. In terms of the real 
and imaginary part of the self energy we can work out 
h(T,uj) to get 



h(T,v) 



(#)' 

, y 



tan 



— — y 
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tan 



2tt 



In 



y 2 + (i-x) 



(13) 



where x = cj— Ei(a;) and y — — £2(^0 with Si (w) [£2 (<*>)] 
the real [imaginary] part of the self energy S(a>). If the 
real part of E is neglected and its imaginary part is as- 
sumed to be a constant T as it would be for impurities in 
an infinite band we can work out the integral in Eq. I|12|) 
and get: 



Note that the second term which deals directly with in- 
teractions between the electrons is of order some scatter- 
ing rate over D/2 while the third which has often been 
emphasized in literature, is of order [fcgT '/ (D /2)] 2 and, 
therefore, can be expected to be much smaller than the 
first. In the oxides, as an example, T is known to be 
of order TM Eq. (O was arrived at independently by 
Benfatto et <dJ£ and by Karakozov and Maksimov. 44 It is 
central to any discussion of single band sum rule. While 
the approach taken in Refs. l43l andfiil are somewhat dif- 
ferent much of the basic content is similar. 

Finally, we consider the superconducting state at T = 
for which in BCS theory 



(15) 



where the gap A k can have d-wave symmetry of the form 
A k = Acos(2(/)). Here A is the gap amplitude and 4> an 
angle along the cylindrical Fermi surface. To keep things 
simple, we start with the s-wave case and obtain 
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(16a) 



Eq. i|16a|) is the difference in KE between superconduct- 
ing (S) and normal (N) state which has increased as ex- 
pected since n k given by Eq. fTK|l populates states above 
the chemical potential while the step function of the nor- 
mal state does not. We have also worked out the differ- 
ence in OS to give 



W 



D 



(16b) 



which has dropped in the superconducting state. Kni- 
gavko et alm& argued that W is the quantity to use when 
discussing the OS. Here we note that W itself shows no 
change with superconducting transition. Taking ratios 
with the normal state in Eqs. (|16afl and (|16b|l we obtain 



W S K E 



2A 
~~D 



In 
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W b - W 
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(17b) 



Note that because In (D /A) is expected to be greater 
than one, the normalized KE change is greater than the 
value of the OS drop in this very simplified model. So 
far we have considered only s-wave. The formulas given 
above hold for d-w&ve with A — > A cos(2</>) and an aver- 
age over angles is required. When this is done 
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and 



W N 



W N _ 1 /2A 



(18b) 



These are useful expressions to understand qualitatively 
the physics but are not believed to be accurate. They do 
represent an extreme case where the OS does not follow 
the temperature dependence of the KE and they differ in 
the superconducting case as well. 



The Nearly Antiferromagnetic Fermi Liquid 
Model 



So far we have not considered interactions yet the ox- 
ides are highly correlated systems. A phenomenological 
approach to correlations is embodied in the Nearly An- 
tiferromagnetic Fermi Liquid model (NAFFL) of Pines 
and coworkers^^Siii and this approach is convenient 
to obtain some information on the effect of interac- 



tions. In this approach the important interactions pro- 
ceed through the exchange of spin fluctuations and the 
imaginary part of the spin susceptibility enters a general- 
ized set of Eliashberg equations. A model susceptibility 
often used is^I 



3m{x(q,w)} 



(u)/w sf )xq 



[1 + C 2 (q-Q)T + (^M?F) 2 



(19) 



The parameters are xq the static susceptibility, Q the 
commensurate antiferromagnetic wave vector (jr/a^/a) 
in the upper right hand quadrant of the first Brillouin 
zone and symmetry related points. £ is the magnetic 
coherence length set at £ = 2.5 a in this paper and 
ujsf is a characteristic spin fluctuation frequency. Fi- 
nally, g is the coupling between charge carriers and 
the spin fluctuations. The Eliashberg equations for 
renormalized Matsubara frequencies w(k, iu> n )i renormal- 
ized energies £(k,iw n ), and pairing energy (gap) in the 
superconducting state 0(k, iu> n ) are^^2^2i5fli£i 



with 



and 



Q(k,iaj n ) = MnirVV X S F (k - k', iw n - iu m ) ^/ lU>m \ , (20a) 
£(Ma,„) = -TVV Agj-(k - k', iu n - iu m ) £k ' +5 k/ ' ' ( 20b ) 
(f)(k,iuj n ) = -TVV A SF (k - k', iu n - iu m ) ' lU)m ^ ; (20c) 



D(k,iuj n ) — uo (k, ioj n ) + [e k + £(k,iw„)] 2 + ^ 2 (k,ia; n ), (20d) 



2 

A^(q,^„) = ? — T 9 n ^\n ,/ v (21) 

i + c (q - Q) + (KIMsf) 

Here q is the momentum transfer and iuj n , iv n are the electron and boson Matsubara frequencies, respectively. Finally, 
the filling (n) is determined from 

(n) = --VV e* + £{k,iu, n ) . (22) 

2 k n>0 w2 ( k ' ,W n) + [ £ k +£(k,iW n )] + 2 (k,ZW„) 



For fixed filling (n) this last equation determines the 
chemical potential which changes with temperature for 
bands which do not have particle-hole-symmetry. 

We solved the Eliashberg equations^ for the two mod- 
els of Table [|] with the single parameter not yet speci- 
fied, g, adjusted to get a critical temperature of 90 K and 
100 K, respectively, for model A and model B. The re- 
sults for W (solid down-triangles) and — Wke/2 (solid 



up-triangles) are shown in Fig. ^ m which they are com- 
pared with the non interacting results. [A 128 x 128 
sampling of the k-space, ak x , ak y £ [0, n] was used which 
accounts for the wiggles. In one case (not shown in the 
figure) the k-space sampling was increased to 512 x 512. 
This made the curves smoother but there were no other 
changes.] In both models the interactions have lowered 
the value of each integral but the temperature depen- 
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FIG. 2: The occupation number rik for selected directions in 
the CUO2 Brillouin zone. Model A of Table[I]was used. Top 
frame: the non interacting case at T — 0. It is included for 
comparison. Center frame: The interacting case at a temper- 
ature T — 20 K. We show normal state (solid line) and super- 
conducting state (dashed line) results. Bottom frame: The 
temperature influence on the normal state for T = 20 K 
(solid line) and T = 150 K (dashed line). 



dence of W and —Wke/2 still largely track each other. 
The role played by the interactions can be traced simply 
in Eqs. QJ and © as changing the probability of occupa- 
tion of the state |k, a) factor n^ a . This factor depends 
on correlations as seen in Fig. [21 where we show results 
for model A. The top frame gives rik (for both spins) 
along r — > X — > M — > T in the first Brillouin zone for 
the non interacting case while the bottom frame shows 
results at T = 20 K (solid curve) and T — 150 K (dot- 
ted curve) in the normal state. We see that the effect 
of interactions is to make the probability of occupation 
of the state |k) non zero even in regions where there is 
no occupation in the non interacting case. Also, there is 
considerable "smearing" of the interacting case distribu- 
tion which changes with temperature and with the onset 
of superconductivity (center frame). On comparing bot- 
tom and center frame we note an increase in temperature 
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FIG. 3: Comparison of normal and superconducting state for 
the optical sum and the kinetic energy. The left hand frame 
applies to Model A of Tableland lusf = 82meV and the 
right hand frame is for the band structure Model B of Table 
and lusf = lOmeV. The dotted line in this frame shows a T 2 
law extrapolation of the normal state data for T > T c to zero 
temperature. Adapted from Ref. IHfl 



of about 100 K in the normal state corresponds roughly 
to the same amount of extra smearing as is due to the 
onset of superconductivity. This smearing in the occu- 
pation factor means an increase in kinetic energy. This 
is expected in a BCS mechanism. In the Cooper pair 
model two electrons are introduced at the Fermi surface 
of a quiescent Fermi sphere. Because of an effective at- 
tractive potential between them they prefer to go into 
a superposition of plane wave states with |k| equal or 
slightly greater than the Fermi momentum, and so 
reduce their potential energy. Although the kinetic en- 
ergy in the process is increased over the free electron 
value of 2e(fci?), the potential is sufficiently reduced to 
compensate for this increase. 

Returning to Fig. ^ several features are to be noted. 
For the case of no interactions the numerical data for 
both W and Wke is well fit by a T 2 law as is also the 
interacting case based on Model A. But this is clearly 
not so for Model B for which the OS and also — Wke/2 
turns up from a T 2 law as T is lowered towards zero with 
turn up onset occurring around T = 100 K. In this case 
the change in W between T — and 200 K is 6.4% which 
is much larger than for the non interacting case (1.5%). 
While the two models A and B have a different Fermi 
surface, dispersion relation, and filling, the most impor- 
tant parameter which gives the deviation from the T 2 
law noted above is the small value of the spin fluctuation 
frequency lusf — lOmeV used in Model B. This implies 
that the thermally activated scattering is larger in this 
model than it is for Model A. Later, when we consider 
simpler versions of our interaction model, we will return 
to this issue where it will become much easier to trace 
these dependencies. 
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In Fig. [3] we compare normal and superconducting 
state for the two models of Table [I] We see that both 
KE and OS (open squares and up triangles, respectively) 
fall below their normal state values (solid squares and up 
triangles, respectively) at the same temperature. This 
corresponds to an increase in KE as the superconducting 
state is entered. The changes are small in all cases. The 
KE in meV has increased by 0.2% for Model A and by 
0.77% in Model B between superconducting and normal 
state at T = 0. For the OS the differences are 0.22% and 
0.6%, respectively, very comparable to what is found for 
the KE but not identical. We note that in the NAFFL 
model once the susceptibility (|2"TJl is specified, supercon- 
ducting solutions with d-wave symmetry result. These 
are not put in by hand as to symmetry or functional 
form. The gap involves a harmonic [cos(ak x ) — cos(ak y )] 
as well as many of the higher harmonics consistent with 
the d x 2_ y 2 irreducible representation of the symmetry 
group for the square Cu0 2 lattice. The solutions are 
far from simple BCS. In d-wave BCS an ansatz on the 
pairing potential Vk,k' is made that it have the form 
?7k^?7k' with 77k ~ [cos(ak x ) — cos(ak y )]. This leads to a 
d-wave gap consisting of the lowest harmonic only. Early 
solutions^ based on Model A with losf — 7.76 meV 
showed that the gap had maximum amplitude at the 
X-point (ir, 0) of the two dimensional Cu02 Brillouin 
zone with A^ax = 33 meV. On the other hand the max- 
imum gap on the Fermi surface was A^fx = 27meV 
at T = 20 K far below its BZ maximum. The value of 
2A™ x /(k B T c ) = 6.27 is quite different from the BCS 
value of 4.28. Also shown on the right hand panel is a 
dotted straight line which indicates the least squares fit 
extrapolation of the normal state data to zero tempera- 
ture from T > T c . It represents a T 2 law. It is important 
to notice, for later reference, that both, the supercon- 
ducting as well as normal state W data are above this 
dotted line for all temperatures T <T C . 



D. The Hubbard Model and Dynamical Mean Field 
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FIG. 4: (Color online) W(Q, T) normalized to its T = value 
as a function of T 2 for the Hubbard model, for Q — 0.05 D 
(top frame) and f2 = 0.75 D (bottom frame). Symbols are 
the results of the DMFT calculations; lines are best fits to 
them. The various symbols give four dopings (x). Adapted 
fromRef.El 



Another approach used to treat strongly correlated sys- 
tems such as the cuprates is Dynamical Mean Field The- 
ory (DMFT). 53 The approach is numerical and based on 
the Hubbard model with hopping t, onsite Coulomb re- 
pulsion U , and chemical potential [i. The half band- 
width D — 8t is taken to be 1.2 eV in what follows with 
U = 3D/2 = 12 1, so that the antiferromagnetic super ex- 
change J = At 2 /U ~ 100 meV. For large values of AU/D 
at half filling the system is a Mott insulator. Away from 
half filling, the weight of the quasiparticle peak at the 
Fermi level denoted by Z and related to the self energy 
E(uj) by 



Z = 



r dsn 




duj 


u=0- 



is non zero but small and is a measure of the metalicity. 



The Hamiltonian used has the form 
H = -t c \<j c 3<y + U Y1 n A n 3i ~ V ( n A +»a). 

(ij)a i i 

(23) 

where a is spin f, J., Ci a (c\ a ) are annihilation (creation) 
operators for fermions of spin a on site i, rii a = c icr Ci a , 
and the sum (ij) is restricted to nearest neighbors only. 
In DMFT the full many body system with strong inter- 
actions is modeled as a single site problem with relevant 
interactions between, say, two electrons at that site plus 
coupling to a bath which represents on average the re- 
maining degrees of freedom. The bath and local problem 
is to be solved for in a self consistent way. The method 
is now well developed and has proved its ability to si- 
multaneously describe low and high energy features in 
Mott systems. Toschi et al£Z have calculated the real 
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FIG. 5: The relative variation of spectral weight between the 
lowest T and 300 K as a function of doping x for various 
cuprates (full symbols) is compared with DMFT calculations 
(open squares) and with the predictions of non interacting 
models [tight binding model (solid line), constant density of 
states approximation (dashed line)] . The dotted line is a guide 
to the eye. The simple inclusion of correlation effects allows 
one to reproduce the observed absolute values with no need of 
fitti ng p arameters. Data for LSCO a re obt ained from Refs.lBil 
and|5l and for BSCCO from Refs. l56l5ll5SkB SCCOi ) and 
from Ref.HI (BSCCO2). Adapted from ReflE3 




FIG. 6: The T = spectral weight W for fi = 0.75D 
as a function of doping x for the Hubbard model (open 
squares) and the experimental observation of Lucarelli et al. 
Ref . |5^ (solid circles) obtained by integrating data on <Ti(w) 
for LSCO. The dotted curve is for reference and gives results 
for U = (no correlations included). Adapted by A. Toschi 
from Ref.E3 



part of the optical conductivity in this model and ob- 
tained the results shown in Fig. 0] for various values 
of the doping x as indicated in the figure which gives 
W{n,T)/W(fl,T = 0) as a function of (2T/D) 2 . The 
left hand frame employs a cutoff f2 in Eq. JIJ of 0.051? 
while the right hand frame is for fl = 0.75D. All curves 



appear to follow a T 2 behavior and the slope of these 
lines becomes smaller with increasing f2. A comparison 
of their DMFT results with experiment is also offered by 
the same authors and we reproduce it in Fig. [3] where we 
show the difference AW = W(T = 0) - W(T = 300 K) 
for fl — 1.5 D renormalized to Wo = W(T = 0) denoted 
by AW/Wo on the figure as a function of doping x (open 
squares). 

Experimental results are for La2- :c Sr ;c Cu04 (LSCO, 
solid circles), Refs. an d l55l and Bi 2 Sr 2 CaCu2 08+x- 
Solid squares are from Refs. l56l57l58l (BSCCOi ) and solid 
triangles from Ref. 1591 (BSCCO2). Also shown for com- 
parison are results for tight binding with no interactions 
(solid line) and in a constant density of states approxima- 
tion (dashed line). It is clear that correlations dominate 
the observed change in W between zero and 300 K to its 
T = value. DMFT provides a reasonable fit to the data 
while tight binding fails badly, giving values that are too 
small. 

In Fig. [S] DMFT results (open squares) as a function 
of doping x are compared with the experimental results 
(sold circles) of Lucarelli et al~ for the optical sum of 
LSCO at T = 0. The theory reproduces well the dop- 
ing dependence of W = W(n = 0.75L>,0). While the 
absolute value of Wq is somewhat underestimated its 
large decrease as the Mott transition is approached is 
well captured by the model calculations. Also shown on 
the same figure (dotted line) are results obtained in the 
limit [7 = (no Hubbard U). We see that in this case the 
observed trend with doping is not reproduced and Wo is 
much too large, particularly at small dopings indicating, 
as expected, that near the Mott transition Wo goes like x 
rather than 1 — x. Finally, we note that the inclusion of 
U (correlation effects) reduces Wo just as we have seen in 
the NAFFL model. We note, however, that this model is 
not well suited to describe doping differences because the 
phenomenological spin susceptibility (jTT?)l really needs to 
be fixed to some experimental observation at each doping 
level, particularly as the Mott transition is approached. 

One thing to note about the DMFT results shown in 
Fig. 0] is that even for the smaller value of the cutoff £1 
on the OS, a T 2 behavior is observed at least for the lim- 
ited data available in the figure. By contrast, we have 
seen m the NAFFL calculations of Figs. Hand CD that de- 
viations from T 2 can occur when the characteristic spin 
fluctuation energy losf is small even when the sum is 
taken over the entire band, i.e.: f2 is large enough to in- 
clude all contributions. Thus, the two calculations differ 
in this important point. We will see later that there is no 
guarantee that a given model for the interactions should 
always give a T 2 law. 

To end we mention related Hubbard model work by 
Maier et ali— They work directly with the KE rather 
than with the conductivity as did Toschi et alJ& and 
use the dynamical cluster approximation (DC A). They 
present results in both normal and superconducting state 
for two doping levels, namely x = 0.05 and x — 0.2. In 
both cases the KE is found to decrease with decreasing 
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temperature even in the superconducting state and this 
decrease, in fact, drives superconductivity. The changes 
are of order a few percent at most. 

As we have been so far mainly interested in the nor- 
mal state we postpone further discussion of this work to 
a later section which deals explicitly with the supercon- 
ducting state. We will also describe the recent work of 
Haule and Kotliar 61 which deals with the t-J model (see 

sec. nvi . 



III. THE EFFECT OF INTERACTIONS IN 
ISOTROPIC BOSON EXCHANGE MODELS 

A. Self Energy in Finite Bands 

Recent studies of finite band effects have revealed that 
the self energy due to impurities or to interaction with 
bosons is profoundly changed by the application of a cut- 
off ±D/2 in a half filled band. 62 i 63 ' fi4 i 65 i 66 i 6 7i 68 In a finite 
band, the normal state self energy due to the electron- 
phonon (or some other boson) interaction is given by 



£(z) = r)^A(z- iuj m )r)(iuj m ) 



+ J dva 2 F(z){[f(T,v- z) + n{T,v)} 
o 

x v (z -v) + [f(T, u + z) + n(T, v)]r)(z + v)} , 

(24a) 

in the mixed representation of Marsiglio et alM> Here, 



\{z) = / dva 2 F(v)- 



2v 



(24b) 



and 



rj(z) = I de 



N {e) 1 
JV o (0) «-€-£(*)' 



(24c) 



These equations need to be solved self consistently as 
£(z) itself depends on itself through Eq. (|24c|) . Here 
z is a complex variable, n(T, v) is the Bose ther- 
mal factor, a 2 F(z) the electron-phonon spectral den- 
sity and N (e)/N o (0) the normalized bare electronic den- 
sity of states taken here to be constant and confined to 
[— D/2, D/2]. This provides a band cutoff in r)(z) of 
Eq. I|24c(l and makes Eq. I|24a(l self consistent. These 
equations are also given by Karakozov and Maksimov 44 
where they are written either in pure Matsubara nota- 
tion or fully on the real frequency axis rather than in the 
mixed notation of Eqs. I|24|) which includes both versions, 
Matsubara and real axis. We have found the mixed no- 
tation more convenient for numerical work. For coupling 
to a single Einstein boson mode a 2 F(uj) = a5(uj — He) 



where He and a are taken in units of D/2. The electron 
mass renormalization is given by A = 2a/ He which is 
dimensionless. For an infinite band at T = we recover 
from Eqs. (|24|l the very familiar result 



£i(u;) = aln 



- n, 



and 



£2(0;) = —ita9 (\uj\ — He) . 



(25) 



(26) 



where 9(x) is the theta function equal to one for x > 
and zero otherwise. For positive values of uj the real 
part £i(w) < and approaches zero as uj 3> He- The 
imaginary part £a(w) is zero till \uj\ = He at which point 
it jumps to a value of ira and then stays constant. Finite 
band effects change this behavior radically. 

In Fig. [7] we show results of Knigavko and CarbotteS& 
for minus the imaginary part [—£2^)] of the self energy 
in frame (a) and its real part [£1(0;)] in frame (b) for 
coupling to a single Einstein phonon at He = 0.1 with 
mass enhancement parameter A = 2 for different temper- 
atures, all in units of D/2. The temperatures are as given 
in the figure caption. The sharper curves correspond to 
the lower temperatures. For T = 0, £i(w) has a loga- 
rithmic like singularity at uj = He as in the infinite band 
case of Eq. (|25|) but now, rather than remain negative 
as it goes towards zero for uj ^> He, £i(w) crosses the 
w-axis and takes on large positive values before dropping 
towards zero from above. Equally different from the infi- 
nite band case, —£2(0;) is not a constant equal to ira [see 
Eq. above u > He but rather drops as we approach 
the bare band edge at uj = 1.0 after which it becomes 
small for uj ~ 1.5 which is where the renormalized density 
of states N(uj) also becomes small. For further discussion 
of finite band features we refer the reader to Refs. |62| to 
l68l as well as to Ref. where a somewhat different den- 
sity of states model is used but this does not change the 
qualitative behavior of the self energy seen in Fig.[7| An 
advantage of Karakozov and Maksimov's choice of DOS, 
however, is that they can get a simple analytic expression 
for the self energy which they also show to be accurate 
by comparing with numerical results based on the full 
Eliashberg equations. 

We return next to the OS of Eq. I|10|) . One can get in- 
sight into the relative importance of the temperature de- 
pendence carried by the thermal factor f(T,uj) in W(T) 
and the temperature dependence solely due to the inter- 
actions in A(e,uj). In Fig. [S] we show results obtained 
by Knigavko et almS, for n(e), Eq. JSJ. Frame (a) shows 
n(e) vs the renormalized energy 2e/D for four values of 
reduced temperatures t = 2T/D, namely 0.0 (solid line), 
0.01 (dashed line), 0.03 (dotted line), and 0.05 (dash- 
dotted line) when only the temperature dependence in 
A(e,uj) is included, i.e.: the f(T,uj) factor in 71(e) was 
excluded. Here we have considered coupling to a sin- 
gle phonon of energy Qe = 0.04 with mass enhancement 
A = 1. It is clear that considerable temperature depen- 
dence comes from this source alone and will influence the 
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FIG. 7: Evolution of the self energy vs frequency dependence 
with temperature for the Einstein spectrum with a = 0.1 and 
Qe = 0.1 (A = 2). Frame (a) is for minus the imaginary 
part, — £2(1^), while frame (b) is for the real part, £i(w). In 
each frame the different curves correspond to temperatures 
t = 2T/D = 0.1 (dash-dotted line), 0.05 (dotted line), 0.01 
(dashed line), 0.001. All energies are in units of D/2. Adapted 
from Ref . |M 



temperature dependence of the OS. In frame (b) of Fig. 
we show results (solid curves) when both sources of T de- 
pendence in n(e) are included and compare with the two 
equivalent curves of frame (a) for t = 0.01 (dashed line) 
and 0.05 (dash-dotted line). It is clear from these fig- 
ures that the temperature dependence of the self energy 
is always important in determining the probability of oc- 
cupation, n(e), of a state of energy e and hence the OS's 
T dependence. 



B. Results for a delta- function boson model 

Results for the OS normalized by a factor of h 2 /(2mb) 
for the tight binding case (|10f> and by h 2 n/m for the 
quadratic dispersion, both denoted by S(t), are given in 
Fig. El reproduced from Knigavko et alm& The parame- 
ters are as shown on the figure and given in the caption. 
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FIG. 8: The probability of occupation n(e) of a state of energy 
e vs normalized energy 2e/D. (a) The results for the case 
when only the temperature dependence of the self energy is 
included. The normalized temperatures are t = 2T/D — 0.0 
(solid line), 0.01 (dashed line), 0.03 (dotted line), and 0.05 
(dash-dotted line), (b) Comparison of the results when the 
complete temperature dependence of n(e) [see Eq. 1101 1 is 
accounted for (solid lines) with the case presented in frame 
(a). The normalized temperatures are t = 0.01 (dashed line) 
and 0.05 (dash-dotted line). (All energies are in units of D/2.) 
A constant density of states model with cutoff ±15/2 is used. 
Adapted from Ref. 03 



The top frame is for a = 0.02 and the bottom frame 
for a = 0.1. The former corresponds to a conventional 
metal while the latter may be more characteristic of the 
cuprates. Two representative values of the boson energy 
£Ie have been chosen to get A = 0.8 (squares) and 1.0 
(triangles) in the top frame and A = 2 (squares) and 
A = 4 (triangles) in the bottom frame. Both tight bind- 
ing (open symbols) and quadratic (solid symbols) dis- 
persion relations are considered. It is quite clear from 
these results that the temperature dependence of the OS 
needs not be quadratic and depends on the size of the 
coupling to the phonons (A value) and also on the dis- 
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FIG. 9: The variation of the optical sum S vs reduced tem- 
perature t = 2T/D for the interaction strengths (a) a = 0.02 
and (b) a = 0.1. The results are for both the quadratic (solid 
symbols) and tight binding (open symbols) bands, as indi- 
cated. In the top frame the mass enhancement factor A = 0.8 
(squares) and A = 1 (triangles), while in the bottom frame 
A = 2 (squares) and A = 4 (triangles). The inset in the 
top frame illustrates the behavior of the optical sum during a 
sudden "undressing transition" at t un dress = 0.02 with 20% 
hardening of the normalized boson frequency Q. and a cor- 
responding reduction of the mass enhancement factor from 
A = 1.0 to 0.8. Adapted from Ref. liul 



persion relation used to describe the electron dynamics. 
An important observation for what will follow is con- 
tained in the inset of the top frame of Fig. |5J The heavy 
solid line follows the temperature evolution of the OS 
when a sharp "undressing" transition is assumed to take 
place at t = 0.02 where a 20% hardening of the phonon 
spectrum occurs. This corresponds to a shift in the mass 
enhancement parameter from A = 1 to 0.8. Such a hard- 
ening of the phonon energy leads to an increase in the 
OS corresponding to a decrease in the KE of the charge 

r, n rr-^rc 40 



C. Results for an MMP spin fluctuation model 

The delta function spectra used in the work of Kni- 
gavko et al. 40 while giving us insight into the mecha- 
nism that leads to temperature dependences in the OS 
is not realistic for the cuprates. If one wishes to re- 
main within the framework of a boson exchange mecha- 
nism one should at the very least choose a different form 
for a 2 F(cu). For spin fluctuation exchange a much used 
choice is the MMP model of Millis et a ;, 47 i 70 i 71 In its sim- 
plest form a 2 F(uj) which will now be denoted I 2 x(to) is 
taken as 



I\H = I 2 



Lo/lo SF 



1 + (lo/lo S f) 2 ' 



(27) 



with / a coupling between charge carriers and spin sus- 
ceptibility and losf a spin fluctuation frequency. Here, 
the coupling I is adjusted to get a superconducting tran- 
sition temperature T c of 100 K when used in the corre- 
sponding Eliashberg equations generalized for a d-wave 
gap symmetry. In this section we concentrate on the 
normal state only. Thus, Eqs. I|24|l apply with a 2 F(ui) 
replaced by Eq. I|27|) . Results for a case with u>sf = 
20meV, I 2 = 0.82, and D/2 = 400 meV are shown in 
Fig. 1101 which has three frames. In the top frame we show 
— Y,2(T, lo) for four temperatures as labeled, in the center 
frame £i(T, lo), and in the bottom frame h(T, —lo) vs lo. 
It is clear that both real and imaginary part of T,(T,lo) 
have important T dependencies which get reflected in 
h(T,ui) and consequently in the OS given by Eq. H10|) 
and more approximately by Eq. l|12f) with the last term 
dropped. In this instance it is simply the area under 
h(T,u>) for negative lo that determines the OS while in 
Eq. (|10fl it is its overlap with the Fermi distribution and 
this has an additional temperature dependence due to 
f(T,u) but it is small. Results for 



S(T) = 2W{T)^ 



(28) 



are shown in Fig. which has two frames. In the left 
hand frame we show results for S(T) vs T 2 up to 40000 K 2 
for our MMP model with lo sf = 20meV and I 2 = 0.82 
while in the right hand frame losf is increased to 82 meV 
with I 2 = 0.655. In both cases the full band width 
D ranges over 2500 (dashed line), 800 (solid line), and 
500 meV (dash-dotted line). It is clear that decreasing 
D reduces the value of S(T = 0) and gives a stronger T 
dependence which, nevertheless, remains near a T 2 be- 
havior for the right hand frame but not in the left hand 
frame. These results are in qualitative agreement with 
those shown in the previous section, Figs. ^ an d El ob- 
tained in a tight binding model with a k-point sampling 
method of the Brillouin zone and the full momentum de- 
pendent Eliashberg equations (|20() . Also shown in both 
frames as a dotted curve for the case D = 800 meV are 
straight lines which allow us to judge better the devia- 
tions from a T 2 law at small T in the curves of the left 
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FIG. 10: Top frame: minus the imaginary part of the car- 
rier self energy — T,2(uj,T) vs lu at the various temperatures 
shown. Middle frame: minus the real part of the carrier self 
energy £i(u;,T) vs uj. Bottom frame: the function h(T,—oj) 
of Eq. 111^ vs uj. The calculation is for an MMP model, 
Eq. with u>sf = 20meV, I 2 = 0.82, and D = 800 meV. 
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FIG. 11: The reduced optical sum S(T) according to Eq. 
vs temperature (T) squared. The left hand frame is for an 
MMP model with lu S f = 20meV and I 2 = 0.82 while 
the right hand frame is for losf = 82meV and I 2 — 0.655. 



exact T 2 law. The actual T dependence of S(T) comes 
from the T dependence of the area under the curves for 
h(T,—uj) (bottom frame of Fig. I1U|) . We see that these 
curves change most at small w as T is increased but there 
are also important changes at large uj and we are sam- 
pling an average of these changes. Nevertheless, the ui 
dependence of a 2 F(uj) or I 2 x{u) at low energies ui will 
determine most strongly the temperature dependence of 
h(T, —uj) at those energies and, thus, plays an important 
role in the temperature dependence of S(T). However, 
even at higher energies there is still significant T depen- 
dence in h(T, —uj). 



D. The non T 



temperature dependence of the 
optical sum 



hand frame. For the right hand frame there is little de- 
viation. Note that both these behaviors are completely 
consistent with numerical results of Fig. [3] and show that 
the upward deviation at small T of the solid curve as 
compared with the dotted (T 2 ) is a robust result not 
dependent on the model used, provided that the spin 
fluctuation frequency is small. Note also that Eq. i|10|) 
includes T variations in h(T, uj) as well as the thermal fac- 
tor f(T,uj), but this latter factor can be neglected. The 
differences are very small confirming our previous claim 
that the temperature dependence in S(T) is mainly due 
to variations in the interaction term at least for these 
cases. Finally, we stress again that the T 2 law found in 
the DMFT calculations for the Hubbard model and of- 
ten confirmed in experiments also arises in the NAFFL 
model with an MMP interaction provided the spin fluc- 
tuation energy is fairly large and D not too small. In 
principle, however, there is no fundamental reason for an 



Benfatto et a/4i have reconsidered the effect of 
interaction on the temperature dependence of S(T) with 
an aim at providing simple analytic expressions that 
would help trace more directly the dependence on the 
interactions. In Fig. 1121 we reproduce some of their re- 
sults for the temperature dependence of the negative of 
the KE, ~{K) = K(T), in meV. The two dotted lines 
at the top of the graph represent the noninteracting case 
(A = 0) and are for reference. As expected, there is little 
difference between K(T) and its zero temperature value 
K(T = 0). The three solid curves are for a system with 
coupling to a boson of frequency uje — 2 meV and A = 1. 
We see that K(T = 0) and K(T = 0) + if Somm (T), the 
Sommerfcld contribution, are not significantly different. 
Thus, Ks omm (T) is not the important contribution to the 
temperature dependence of the complete KE [heavy solid 
line, denoted K(T)\ which develops a linear temperature 
dependence for T ~ 20 K. The same holds for the case 
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FIG. 12: (Color online) The negative of the KE, —(K), vs 
temperature showing the breakdown for the various contri- 
butions. The solid curves are for a small boson frequency 
uie = 2meV with A = 1. The two thin lines show how the 
zero temperature value is modified by the Sommerfeld term 
by a very small amount. A similar size contribution (actually 
a bit larger) is illustrated for the noninteracting case by dot- 
ted curves, where only the Sommerfeld term is responsible for 
the temperature variation. Finally, for higher boson frequen- 
cies, the dashed curves illustrate the amount of temperature 
variation due to the Sommerfeld term compared with the rest. 
The band width W = D = 1 eV. Adapted from Ref. |H 
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FIG. 13: (Color online) The partial spectral weight integrated 
up to various frequencies as a function of temperature. Be- 
low 16 000 cm" 1 there is an increase in spectral weight as the 
temperature is lowered signaling a line narrowing on this fre- 
quency scale. Below T c there is strong loss of spectral weight 
to the superconducting condensate. There is no evidence of 
any precursors to superconductivity at 67 K. In the inset we 
show N e ff(u) at 295 K. Adapted from Ref. 171 



lue = 25meV (three dashed curves) but the variation in 
T of K(T) in this case is less pronounced and closer to a 
T 2 law as was found for the MMP spectrum. The band 
width in all three cases was W = D = 1 cV. Benfatto et 
alJl trace these T dependences in detail using analytic as 
well as numerical techniques. They concluded that both 
the real and imaginary parts of X contribute. Karako- 
zov and Maksimo\i^ also find deviations from a T 2 law 
and provide an approximate analytic formula for these 
deviations in a particular case. For different models of 
a 2 F{uj) other T dependencies are possible. This means 
that in principle one can learn about some features of 
this underlying a 2 F(uj) from a study of the temperature 
dependence of the OS, but the correspondence is not nec- 
essarily simple or unique. On the other hand, a complete 
study of the temperature and frequency dependence of 
the optical self energjiZ£ as is now done routinely, pro- 
vides much more detailed information than does the OS. 
Recall that h(T,u>) is dependent on an average over the 
real and imaginary part of S, the quasiparticle self en- 
ergy. This information is accessible directly from angular 
resolved photoemission spectroscopy (ARPES) at each 



point in the Brillouin zone separately. 



While many experiments give a T 2 dependence in the 
normal state within the precision of measurement there 
is some evidence for other dependences in the cuprates. 
In Fig. EH we reproduce data from the paper of Hwang et 
a\32k for the partial sum to uj c in units of carriers per Cu 
atom denoted by N e ff(T). What is shown is the differ- 
ence N e f f (T) — N e f f (300 K) as a function of temperature 
T to 300 K. The data is for underdoped YBa 2 Cu306.5o 
(YBCO6.50) Ortho II material of high quality and pu- 
rity. There are 5 values of u c , namely 1000, 5000, 10 000, 
16 000, and 20 000 cm -1 . It is clear that the tempera- 
ture dependence of N e f f (T) is not necessarily quadratic 
in T for the partial sums and also for the case most rel- 
evant for our discussions with tu c = 16 000 cm -1 . These 
samples are, however, underdoped and a pseudogap is in- 
volved which could have some effect on the temperature 
dependence of the OS as we will discuss later within a 
phase fluctuation model for the pseudogap and alterna- 
tively a D-density wave model. 
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FIG. 15: The charge carrier-spin spectral density I 2 x{u) de- 
termined from optical scattering data of YBCO6.95 at vari- 
ous temperatures. Solid gray curve T — 90 K, dash-dotted 
T = 80 K, dashed T = 60 K, dotted T = 40 K, and black solid 
T — 10 K. Note the growth in strength of the 41 meV opti- 
cal resonance as the temperature is lowered. Adapted from 
Ref.El 
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FIG. 14: Experimental values of the 06-plane spectral func- 
tion pl defined in Eq. H29H . Adapted from Ref. I34 



IV. THE OPTICAL SUM IN THE 
SUPERCONDUCTING STATE 

A. Superconducting optical sum including of 
inelastic scattering collapse 

Next we return to a more detailed look at the super- 
conducting state. In Fig. El we reproduce the original 
results of Molegraaf et a/. 59 for two samples of BSCCO 
as they were presented by van der Marel et idml The top 
frame is for optimally doped and the bottom frame for 
an underdoped sample. What is shown is 

2 n 

PL = C^ [ duUe[a(cj)] (29) 
Tre J 

—Q 

in meV with C a constant to be specified shortly. It 
is striking that in these two samples the OS appears to 
increase faster in the superconducting state than in the 
underlying normal state extrapolated to low tempera- 
tures (dotted line). This was interpreted by Molegraaf et 
alm^ as indicating a decrease in KE in the superconduct- 
ing state as compared with the underlying normal state 
which also shows a decrease in KE as T — > but by less 



than in the superconducting state. This behavior clearly 
goes beyond ordinary BCS theory and also Eliashberg 
theory as formulated for phonons. 

There remains some controversy about the interpreta- 
tion of such data with some authors finding the opposite 
effect 2& We will not go into these details here but re- 
fer the reader to some relevant literature£iZ4i££ and note 
that Santander-Syro^ confirm the basic results of Mole- 
graaf et al~ See also the recent work of Carbone et alm& 

We have already seen in Fig. [5] [inset in frame (a)] 
for the normal state, that undressing effects associated 
with a hardening of the phonon spectra at a particular 
onset temperature can give an OS curve that behaves 
very much like those seen in Fig. 1141 While in princi- 
ple phonon frequencies can shift as a result of the onset 
of superconductivity these effects are small and usually 
negligible. For an electronic mechanism, however, we 
have already discussed in the introduction the idea of the 
collapse of the inelastic scattering as superconductivity 
sets in which leads directly to a peak in the real part of 
the microwave conductivity at some intermediate tem- 
perature below T c . Within a spin fluctuation mechanism 
this translates into a hardening of the spin fluctuation 
spectrum and the modification of I 2 x{u) of Eq. 1|27H as 
the d-wave superconducting gap opens. An analysis of 
the optical scattering rates in YBCO6.95 by Carbotte et 
al& showed that there is a reduction in I 2 x{oj) at small 
lu as the temperature is reduced. Results obtained by 
Schachinger et alS^ for the temperature evolution of 
7 2 x(cj) in YBCO6.95 are reproduced in Fig. ^1 We see 
a reduction in spectral weight at low u> as well as the 
growth of an optical resonance at 41 meV, the energy 
of the spin one resonance seen in neutron scatterings^ 
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FIG. 16: Temperature dependence of the conductivity <Ti(w) 
at microwave frequency u) = 0.144 me V in YBCOe.gs. The 
open circles are results based on our model spectral density 
(see Fig. I15H obtained from inversion of optical conductivity 
data and the dashed line is based on an MMP model with 
low frequency cutoff. The solid triangles include impurities 
with the solid line from Ref. The solid squares represent 
experimental data by Bonn et aZp— Adapted from Ref. WH 



When these spectra are used in a generalization of the or- 
dinary Eliashberg equations that includes the possibility 
of d-wave symmetry with gap A(</>) = Ao cos(20) where 
<t> is an angle on the Fermi surface taken to be cylindri- 
cal, Schachinger and Carbottefi 8 *^! find excellent agree- 
ment for the microwave data of Hosseini et alm^ taken at 
five different microwave frequencies on high quality sam- 
ples of YBC0 6 . 99 . In Fig. [Tgl reproduced from Ref. |hI 
we show older results obtained by Bonn et alS^ fit by 
Schachinger et al. 83 by simply applying a low frequency 
cutoff to an MMP formSSi The data are given as the 
solid squares while the theoretical calculations without 
impurities (open circles and dashed line) and with impu- 
rities (solid triangles and solid line) are shown for com- 
parison. Here t + gives the impurity scattering rate in 
Born approximation. Open circles and solid triangles are 
based on the spectra of Fig.ll5lwhich include the 41 meV 
spin resonance. The dashed and solid lines are from ear- 
lier calculations based on an MMP spectrum with appli- 
cation of a low frequency cutoff without resonancei 83 i 84 
We see that such a procedure gives results that are very 
close to those based on the more exact results of Fig. \U>\ 
Calculations in the model of Eqs. (|20|l of the OS and 
KE of optimally doped BSCCO have been carried out by 
Schachinger and Carbotte^S who simulated the expected 
hardening of the spin fluctuation spectra in the super- 
conducting state by applying a low frequency cutoff to 
the interaction of Eq. 1)21(1 . This cutoff varies with tem- 
perature. It is zero at T c and has a maximum amplitude 
of 23 meV at T = with intensity changed according to a 
mean field BCS d-wetve order parameter temperature de- 
pendence. Their results are shown in Fig. 1171 The quan- 
tity pl according to Eq. I(29|l on the vertical axis is in meV 
and the horizontal scale is T 2 in K 2 . The heavy solid line 
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FIG. 17: The optical sum of optimally doped BSCCO as a 
function of the square of the temperature for the band struc- 
ture Model A of Table U with different values of ujsf- Also in 
one case a low frequency cutoff is applied to the spin suscep- 
tibility. Note the significance of the T 2 -variation on the value 
of ujsf- The thick solid line represents experimental normal 
state data of Molegraaf et al*— also shown in the top frame 
of Fig. E| Adapted from Ref. ET1 



are normal state experimental data of Molegraaf et al.*^ 
also shown in the top frame of Fig. ll4l of this review. The 
theoretical results for the OS are all based on Model A of 
Table[I]but were scaled upward by a factor [C of Eq. J2jJ] 
of approximately two to fit the data. The spin fluctua- 
tion frequency ujsf was also adjusted to improve the fit; 
ujsf = 13meV is best. To reduce the discrepancy in 
absolute value of the OS the value of the nearest neigh- 
bor hopping parameter t in our band structure model 
would need to be increased. This would, however, also 
decrease the sensitivity of the OS to temperature varia- 
tions and so ujsf would have to be adjusted downwards 
as well. Markiewicz et alr& have suggested significantly 
larger values of nearest neighbor hopping t for the bare 
band structure of BSCCO than used to fit ARPES data. 
As is well known, there can be a factor of two or more. 
We have used the dispersion relation of Table I of Ref. Is^l 
(specifically we used t = 360 meV, t' = — 100 meV, and 
(ri) = 0.4. We also included further neighbors, namely 
t" = 35meV and t'" = 20meV) and found a value of 
W(T = 0) of 230 meV well above the experimental value 
with ujsf = 8 meV giving the proper temperature depen- 
dence. It is clear that a value of t intermediate between 
that of Markiewicz et al$^ and the one in Table [Q for 
Model A would be needed to get agreement with the OS 
data without any adjustment and a ujsf value between 8 
and 13 meV. Schachinger and Carbotte5 Q did not attempt 
such a fit as their aim was not to fit any particular case 
but rather see how a low frequency cutoff applied to a 
spin fluctuation model changes the OS. The solid squares 
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in Fig. El show the normal state results while the solid 
triangles are in the superconducting state. These are ob- 
tained without cutoff and so the KE has increased with 
respect to its normal state. The open squares and tri- 
angles show results for the same two cases but now the 
low frequency cutoff is applied to simulate the harden- 
ing of the spectrum as T is reduced below T c . We see 
a large decrease in KE as compared to the case without 
cutoff and the superconducting state would then show a 
decrease in KE as compared to the normal state with- 
out cutoff in agreement with the data of Fig. El It is the 
modification of the underlying interactions that have lead 
to this effect. It is clear that a modest hardening of the 
spectra consistent with the changes seen in the electron- 
boson spectral functions of Fig. EJ w °uld be sufficient to 
explain the data in Fig. El Of course, in a complete the- 
ory as yet not attempted, it would be necessary to find 
an interaction which can give not only the correct OS 
but also the microwave peak and all the other properties 
of the superconducting state. Finally, we note that in 
Fig-Elthe OS shows a non T 2 behavior at low tempera- 
tures even without the low frequency cutoff which has the 
effect of making it more pronounced. If we did not have 
the normal state data below T c and simply extrapolated 
the normal state data above T c with a T 2 law to T = 
we would have to conclude that the OS is higher in the 
superconducting state than in the extrapolated "normal 
state" yet in this case the mechanism for pairing is not 
exotic in any way, i.e.: there is no low frequency cut off. 
(See also Fig. [21) With the introduction of a low fre- 
quency cutoff this effect becomes even more pronounced 
and for the case shown the rise in the superconducting 
state is larger than seen in the experiments of Fig. El We 
caution the reader, however, that while this rise is seen 
by several experimental groups and is not controversial, 
its actual size is. 



B. The Temperature Dependent Scattering Time 
Model 



Recently Marsiglic 8 — has given calculations that pro- 
vide support for the results of Fig. II 71 using a related but 
much simpler model. He notes that in an infinite band 
for an elastic scattering rate T the probability of occu- 
pation of the state |k) at finite temperature is given by 
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FIG. 18: (Color online) Minus the kinetic energy, —(T xx ), vs 
temperature squared for various degrees of elastic scattering. 
Below T c we change the elastic scattering rate smoothly to 
zero [as Fq(T/T c ) 4 ] as the temperature is lowered. Adapted 
from Ref. |8u 



function. Note that at zero temperature 



£k 2 _, ( 2Ek 



(32) 



which, in the limit T — > 0, reduces to the well known ex- 
pression (1 — e k /S k )/2, Eq. (JT3J introduced in Sec. Ill Bl 
It is clear from Eq. I|32|l that both, the appearance of 
a gap in i5 k and the elastic scattering T smear n,s(e k ). 
From these expressions for n(e k ) it is easy to work out 
the average KE which Marsiglio denotes by —(T xx ). He 
uses a tight binding band and a constant Tq = 1/r = 0, 
2, 5, 10, and 15meV from top to bottom in Fig. El For 
the superconducting state To is changed to T(T) with 
r(T) = r (T/T c ) 4 which is observed in the work of Hos- 
seini et al£^ This introduces phenomcnologically the idea 
of the collapse of the scattering rate. We see in Fig. El 
that for 1/r = 5 meV and larger the KE decreases below 
its normal state extrapolation as the system becomes su- 
perconducting. The change is of the order of a few per- 
cent for the largest 1/r considered and is large enough 
to explain the measured KE changes in the underdoped 
cuprates. 



njv(e k ) = xsmip - + —— + i- 



f , i £ k i^ . A r E k 

ns et = \smvj — h i 

M ' 2 £ k 7T v \2 4ttT 2ttT 



(30) 



2 7T V 2 47rT 27TT 

for the normal state and, to a good approximation, by 



(31) 



for the superconducting state with £7 k = yi^ + A k , with 
A k the superconducting gap and ip(z) is the digamma 



C. The model of Norman and Pepin 

The mechanism described above which leads to an OS 
increase as superconductivity sets in has some common 
elements with the ideas of Norman and Pepin^LS al- 
though there are also important differences. These au- 
thors use a less fundamental, more phenomenological 
approach based on ARPES^S and optical conductivity 
data from which they construct directly a model for 
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the self energy For the normal state they begin 

with a constant, frequency independent scattering rate 
r = QmE(oj) (taken to be of order 100 meV) to simulate 
the very broad (incoherent) line shapes seen in ARPES at 
the anti nodal point (0, ir). In the superconducting state 
a sharp coherence peak appears which signals increased 
coherence. To model this fact a low frequency cutoff ojq 
is applied to the imaginary part of S(w) making it zero 
for lu < luq and Tq above. The value of luq is chosen 
as the energy of the spectral dip feature seen in ARPES 
in the superconducting state. Here, as in the model de- 
scribed in Sec. IIVBI the underlying normal state on top 
of which superconductivity develops effectively has re- 
duced scattering, i.e.: is more coherent which is a critical 
additional feature not contained in ordinary BCS and 
this leads to a reduction in KE and, consequently, in an 
increase of the OS. To arrive at their final estimates for 
the KE increase Norman and Pepin include further com- 
plications in their model such as the anisotropy observed 
in scattering rates22*2ii22ii£ (r k instead of Tq) going from 
anti nodal to nodal direction as well as an uj dependence 
(proportional to a momentum independent parameter a) 
in the self energy modeled on the Marginal Fermi Liquid 
model 94,95 with parameters determined from optical data 
on scattering rates. We reproduce in Fig. ED their final 
estimates for the KE change —AEk between supercon- 
ducting and normal state. The left hand frame shows the 
optical scattering rate denoted l/r(w) for four BSCCO 
samples from Puchkov et al£& which they use to fit pa- 
rameters. The right hand frame presents the results for 
the change in KE (—AEk) associated with the forma- 
tion of the superconducting state in meV as a function 
of doping x. The solid circles give the theoretical results 
and the open squares and open diamonds experimental 
data from Santander-Syro et al& and Molegraaf et al.^L 
respectively. The theoretical estimates are deemed rea- 
sonable but not accurate. On the overdoped side the OS 
behaves in a conventional fashion while for the under- 
doped side it is anomalous representing a lowering of KE 
as compared with what is expected in conventional BCS. 
Note that the lowest open square in this graph shown as 
being zero for the overdoped sample represents an early 
not very accurate estimate. In Ref.|57]it is actually nega- 
tive rather than zero bringing it closer in agreement with 
the lowest solid circle (theory). This point was further 
discussed by Deutscher et alm^. 



D. Additional data on kinetic energy changes 

We reproduce in Fig. [201 the data of Santander-Syro 
et o/^i for their overdoped sample as presented by 
Deutscher et al^ What is shown is the optical spec- 
tral weight in units 10 6 Q cm~ 2 as a function of T 2 . 
The solid circles are in the normal state and the open 
circles in the superconducting state below T c = 63 K. A 
clear T 2 law is noted above T c and a reduction in spec- 
tral weight below, as expected in ordinary BCS theory. 
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FIG. 19: Top frame: The optical scatte ring rate l/r(u) vs 
ui for various BSCCO samples from Ref. -96 (OD overdoped, 
OPT optimal doped, UD underdoped). Bottom frame: Cal- 
culated sum rule violation (—AEk) vs doping x (solid circles). 
The curve is T c . Also shown in this frame are the experimen- 
tal results (open squares from Ref. 15^ . open diamonds from 
Ref. l59f) . The theoretical doping trend in (b) is due to the 
increasing offset in 1/r seen in (a). Adapted from Ref. WA 



For the optimally doped sample (not shown) the nor- 
mal state temperature dependence is closer to linear than 
quadratic and the underdoped sample shows a very flat 
region before entering the superconducting state. This 
demonstrates once more that there is as yet no strong 
consensus in the literature as to the T dependence of the 
normal state. (Please see also Ref. l98h The same paper 
also analyses the data in terms of KE change for three 
samples UND70K, OPT80K, and OVR63K, as shown in 
Fig. [21 reproduced from Deutscher et al~ The horizon- 
tal axis is p — p pt where p is the charge per Cu atom 
related to T c by 



T, 



c.opt 



= i - 86.2 (p - Popt y 



(33) 



with T C)0pt the maximum critical temperature for the 
Bi2212 seriesiS 9 . What is clear from this figure is that 
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T 2 [10 4 K 2 ] 

FIG. 20: Spectral weight of the overdoped Bi2212 sample, 
integrated up to leV, plotted vs T 2 , from Ref. I57L Closed 
symbols: spectral weight in the normal state, open symbols: 
spectral weight in the superconducting state, including the 
weight of the superfluid. The errors in the relative variations 
of the spectral weight are of the size of the symbols. Adapted 
from Ref. |U 




P-Popt 

FIG. 21: Change AEkin of the KE, in meV per copper site 
vs the charge p per copper site with respect to p op t [Eq. 11331 1 
in BSCCO. Full diamonds: data from Ref. 58, high-frequency 
cutoff 1 eV. Open circles: data from Ref. I5SL high frequency 
cutoff 1.25 eV. Error bars: vertical, uncertainties due to the 
extrapolation of the temperature dependence of the normal 
state spectral weight down to zero temperature; horizon- 
tal, uncertainties resulting from T c /T Ctrnax through Eq. 1331 . 
Deutscher et alw— took T c ^ ma x = (83 ± 2) K for films and 
(91 ± 2) K for crystals. Adapted from Ref. I35L 



there is a smooth crossover from standard behavior on 
the overdoped side to anomalous behavior on the under- 
doped side. The very recent data of Carbone et alM 
lends further support to this conclusion. 



E. Cluster dynamical mean field of the t—J model 

Recently Haule and Kotlias^i calculated the optical 
conductivity of the t-J model within a cluster DMFT 
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FIG. 22: (Color online) The difference between the super- 
conducting and normal state energies as a function of tem- 
perature. The following curves are shown: up-triangles - 
Ekin-s — Ekin-N; squares - E x -s — E x -n; diamonds - 
Etot-s — Etot-N', circles - [is — (in- Adapted from Ref. l6ll 



(CDMFT). (Please see also earlier work by Maier et alM 
based on the Hubbard model.) Haule and Kotliar^i treat 
the temperature and doping dependence and address the 
issue of the change in KE of the holes when superconduc- 
tivity sets in. They find, in agreement with experiment, 
that on the overdoped side the KE makes a negative con- 
tribution to the condensation energy as in conventional 
BCS theory but that there is a crossover to the opposite 
case on the underdoped side. This is shown in Fig. l2~2l 
reproduced from Ref. 0. What is shown is the differ- 
ence between the superconducting and normal state en- 
ergies (E) as a function of temperature, both in units of 
t (nearest neighbor hopping). The up-triangles give the 
KE contribution (Ek% n ), the squares are the super ex- 
change (E x ), and the diamonds the total energy (E to t)- 
We see a change in sign in the KE contribution upon 
condensation as we go from the overdoped to the un- 
derdoped regime. (8 denotes the doping.) These results 
contrast with earlier work of Maier et alm& also based on 
the Hubbard model. There, the KE is found to be lower 
in the superconducting state than in the normal state for 
both values of doping presented, namely 8 = 0.05 and 
8 = 0.2 (overdoped). In the underdoped case the po- 
tential energy increases slightly above its normal state 
value below T = T c . For the overdoped case it does de- 
crease very slightly but plays a much smaller role in the 
condensation than the corresponding drop in KE. 

We end this section by mentioning two related 
workgiogjioi based on the negative U Hubbard model 
which has been used successfully to describe the BCS 
- Bose Einstein (BE) crossover. Toschi et alw^ employ 
DMFT and Kyung et alM± a cellular DMFT and ob- 
tain very similar results. Both normal and superconduct- 
ing states are considered. Both groups find a change of 
sign in the KE difference between superconducting and 
normal state as one goes from underdoping (anomalous) 
to overdoping (conventional BCS behavior) as seen in 
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FIG. 23: (Color online) Low-temperature behavior of the nor- 
malized kinetic energy i?kin/ D (upper panels) and the normal- 
ized potential energy E po t/\U\ (lower panels) as a function of 
the normalized temperature T/D. Here U is the Hubbard 
attraction. The c ritica l temperature is marked by arrows. 
Adapted from Ref. llOOt 



Fig. [2] In Fig. [23 we reproduce from Ref. llOOl results for 
KE (-Ekin) m units of the band width D vs the normalized 
temperature T/D (upper frames) and for the potential 
energy E pot in units of \U\ (lower frames). Three values 
of U are considered, U — 0AD in the left frames and 
U = 1.2 D and 3.21? in the frames on the right. These 
values correspond to BCS, intermediate, and BE (Bose 
- Einstein) regimes, respectively. The frames on the left 
show conventional behavior with a small increase in KE 
in the superconducting state and a decrease in potential 
energy. On the other hand, in the frames on the right the 
KE decreases while the potential energy increases which 
is referred to as KE driven superconductivity. 



V. MODELS OF THE PSEUDOGAP STATE 



A. The Preformed Pair Model 
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FIG. 24: Kinetic energy per bond, {k x ), as a function of tem- 
perature for non interacting tight-binding electrons (TB), the 
BCS solution (BCS), and the phase fluctuation (PP) model 
for jj, — ((n) = 1). The large vertical arrows indicate the 
increase in KE upon pairing, relative to the free tight-binding 
model, and the small arrows indicate the additional increase 
due to phase fluctuations. This additional phase fluctuation 
energy rapidly vanishes near T c = Tkt, which causes the 
significant change in the OS upon entering the superconduct- 
ing state at Tkt = 0.1 1. Note that the thick line follows 
the actual KE encountered in this model, when going from 
the pseu dogap to the superconducting regime. Adapted from 
Ref.lTol 



and a pairing term 



IS i 



(35) 



where 8 connects the site i to its nearest neighbors only 
and (ij), as before, is limited to nearest neighbors as well. 
Furthermore, 



1 



,t J J J 



(36) 



Its average = A exp (i&is) with A the gap ampli- 

tude is assumed to have d- wave symmetry and the phase 



Another model for superconductivity in the oxides is 
the preformed pair model, 102 ! 103 ! 104 ! 105 ! 106 ! 107 The idea 
is that the pairs form at a temperature T* , the pseudo- 
gap formation temperature. The superconducting state 
emerges from the preformed pair state at a lower tem- 
perature T c when phase coherence sets in. The pseudo- 
gap regime is then due to superconducting phase fluctu- 
ations. In a recent paper Eckl et aL— considered the 
effect of phase fluctuations on the OS, i.e.: the KE in the 
pse udoga p regime. [See also the related work of Kopec, 
Ref. Il09t ] They start with a Hamiltonian which contains 
two terms, the KE 

K = ~t ( C *W + 4 c -) (34) 



q s = I +</>i+s)/2 for 8 along x, 
1 {cpi + <f> i+s ) /2 + 7T for 8 along y. 

The phases are assumed to fluctuate according to classi- 
cal XY free energy. In this model the Kosterlitz-Thouless 
transition Tkt is identified with the superconducting 
transition temperature T c and the mean field tempera- 
ture Tmf at which the pairs form, is identified with the 
pseudogap temperature T*. They take Tkt — T mf /4 
with Tkt = O.K. In Fig. 122 reproduced from Ref.lTol 
we show the KE per bond, (k x ) as a function of the re- 
duced temperature T/Tkt- Here 

(k x ) = -tJ2 (4ci+*,* + 4+*.*^) ■ (38) 



21 



0.004 




T/t 



FIG. 25: Kinetic energy contribution from phase fluctuations 
S(k x ) = {k x )pp — (k x )BCS- One can clearly see the sharp 
decrease of KE near the Kosterlitz-Thouless transition Tkt = 
0.1 1 = T c . AEk gives an es tima te of the kinetic condensation 
energy. Adapted from Ref. Il08t 



The dashed curve is the result of the simple tight-binding 
Hamiltonian [K of Eq. I|34[l only] for zero chemical po- 
tential and (n) = 1. The light solid line is the result of 
BCS mean field and the heavy solid line the result ob- 
tained by taking into account phase fluctuations within 
the preformed pair model. The mean field BCS conden- 
sation increases the KE above its tight-binding value and 
the phase fluctuations provide an additional KE increase 
which vanishes at T c — Tkt and this causes significant 
change in the OS as the superconducting state is entered 
at Tkt = 0.1 1 in this model calculations. The KE gain 
from the phase fluctuations S(k x ) = (k x ) PP - (k 

x) BCS 

is shown in Fig. 1251 which shows a sharp decrease in KE 
near the Kosterlitz-Thouless transition and AEk gives an 
estimate of the KE condensation oc 0.003 1 of the same 
order as was found in other models. However, this is 
not the KE change between normal and superconducting 
state at T = as discussed previously. Rather it is a 
change in KE due to the suppression of phase fluctua- 
tions present in the pseudogap state. At T c = Tkt the 
phases become locked in and 6(k x ) consequently becomes 
zero. The change in KE on formation of the Cooper 
pairs is never explicitly considered in this model as the 
pairs are assumed to form at a much higher temperature 
T = T* = OAt. 



B. The D-density Waves, Competing Order Model 

There are other very different models that have 
been proposed for the pseudogap state which exists 
at temperatures above the superconducting state in 
the underdoped regime. One proposal is D-density 
wavea 110 ! 111 ! 112 ! 113 . 114 . 115 ! 116 ! 117 ! 118 ! 119 (DDW) which 
falls into the general category of competing interactions. 
The view is that a new phase, not superconducting, 
but having a gap with d-wave symmetry, forms at 



T* and the superconductivity which arises only at 
some lower temperature T c is based on this new ground 

g ^ a ^ e | 110 4 lll 4 112 4 113 4 114 i 115 4 n6 4 m^ 

The DDW state breaks time reversal symmetry because 
it introduces bond currentsiS2ii2S with associated small 
magnetic moments and a gap forms at the antiferromag- 
netic Brillouin zone boundary. The mean field DDW 
Hamiltonian is 



H 



(39) 



where the sum on k ranges over the entire Brillouin zone 
and Q is the commensurate wave vector (ir/a, ir/a). Here 
D k = Dq [cos^a) — cos(k y a)] /2 with Do the gap ampli- 
tude. Within a mean field approximation the probability 
of occupation of the state |ker) is given bj»i2£ 



nka (T) 



2D k 



{D k [/(T,£ +k ) + /(T,£- k )] 



-e k [/(T,£ +k )-/(T,£_ k )]}, 



(40) 



where £± k = — /j,±D k and D k = \J t\ + D 2 ,. The number 
density is 

^ = 1 E [/(T,£-k) + /(T,£ +k )], (41) 

keMBZ 

where the sum over k is restricted to half the original Bril- 
louin zone, i.e.: to th e magnet ic Bri lloui n zone (MBZ). 
The KE from Refs. EM [lH and [l27| is denoted by 
W(D, T) and is 



w{d,t) = - ^M/me+k) 

keMBZ k 



/(T,e_ k )]. (42) 



Here and later in this section we give formulas only for the 
simpler case of t' = 0, i.e.: only nearest neighbors. The 
results to be presented, however, are based on straightfor- 
ward generalizations that include second nearest neigh- 
bors hopping. 

Aristov and Zeyhei— I calculated the optical conduc- 
tivity in the DDW model with and without vertex correc- 
tions to the usual current operator. The previous work of 
Valenzuela et alM^ was without vertex corrections as is 
the more recent work of Gerami and Nayak—- who, how- 
ever, consider the effect of anisotropic scattering. These 
authors use a current operator which includes a term pro- 
portional to the gap velocity V k D k which is introduced 
so as to ensure charge conservation as is discussed in more 
detail by Benfatto et al. 126 Here we follow Aristov and 
Zeyhe»i22i and show in Fig. their results for the optical 
spectral weight as a function of temperature T in units of 
the hopping parameter t. The solid curve is with and the 
dash-dotted curve is without vertex corrections. We see 
that while vertex corrections have changed the magnitude 
of the OS they have changed much less its temperature 
dependence which is somewhat more pronounced in the 
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vertex corrections. In Fig. [57| we reproduce results of 
Aristov and Zeyheni^S for the conductivity cr(u>) vs fre- 
quency w for t = 0.25 eV, t' = 0.076, T = 0.001, n = 0.4 
and a scattering rate of 1/t = 0.01. It is seen that both 
the Drude region (intra-band) and the inter-band tran- 
sitions which set in at higher lu > 0.2 are larger when 
vertex corrections are included but that beyond u) > 0.3 
the two curves merge and the vertex corrections are no 
longer important. 

In the continuum limit of the DDW model Valenzuela 
et alw2& have obtained useful analytic results for intra 
(Drude) and inter (verticle transitions) separately. Ver- 
tex corrections are neglected and the limit of zero impu- 
rity scattering was taken. They find 



&intra 



-5(ui)ire 2 



j dv 




dv 





giiy), (43) 



FIG. 26: Restricted optical sums with (solid line) and without 
(dot-dashed line) vertex corrections, and Drude weights with 
(dashed line) and with out (dotted line) vertex corrections. 
Adapted from Ref. Il3cl 



solid curve. In both cases the OS decreases with decreas- 
ing temperature as was the case in BCS theory. Also 
shown in Fig. |23 are results for the Drude weight sepa- 
rately, dashed and dotted lines with and without vertex 
corrections. We see that it is also strongly enhanced by 



with 



9l (u) = (hv F ) 2 N(0)- 



tt/2 



n 2 

cos 2 0, 



(44) 

which can be written in terms of elliptic integrals. Here, 
vp is the Fermi velocity. Note also that the Fermi factor 
f(T, v) has [i — and the chemical potential has been 
transfered to g\(y). For the inter-band 



(j inter (T,uj) = (hvp) N(0) — 

7T 



1 f2D, 



tt/2 



d6 



1 - 



'2A> 



B(T,cj) 



(45) 



r 



where B(T,lu) is a universal thermal factor 



B(T,u) 



sinh(/3w/2) 



to cosh(/3/ii) + cosh(/3w/2) 



= 



1 



ksT 1 
(46) 

which at T = becomes proportional to a theta func- 
tion Q(uj — 2|jLi|). This provides a low energy cutoff to the 
inter-band transitions. When impurity scattering is in- 
cluded the delta function 5(cu) in Eq. I|43|) broadens into 
a Lorentzian and temperature leads to the overlap of the 
two contributions as seen in Fig. 1271 

We note that the kinetic energy shows the same trend 
in its temperature dependence as does the OS of Fig. 1261 
This is shown in Fig. 1281 which we took from the work of 
Benfatto et alm^ and where it is denoted by W(D, T) as 
the dash-dotted curve for the parameters shown in the 
figure and described in the caption. We see the same 
decreasing trend with decreasing temperatures as for the 
OS. There are two other curves, the dashed one is for a 



pure tight-binding band with no DDW and is included for 
comparison. The solid curve is for W DDW (D, T) given 
by 



MBZ 



w 



DDW 



£ k [/(T,£+k)-/(T,£_ k )] 
(47) 

(again in units of ire 2 /h ) in the notation of Benfatto 
et alw2&J2L and requires explanation. It was obtained 
from a current operator which was modified to ensure 
charge conservation without the need for vertex correc- 
tions. This current operator has been used in other 
works^^i as well. As noted by Aristov and Zeyhei— I, 
however, such a procedure tends to overestimate the con- 
ductivity at higher frequencies and the OS now shows an 
increase as T is decreased. 

We can also add a mean field BCS term 



H' 



-kJ.CkT 



h.c. 



(48) 
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FIG. 27: Optical conductivity cr(ui) with (solid line) and with- 
out (dashed line) vertex corrections. Adapted from Ref. Il3u1 . 
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FIG. 28: (Color online) Optical spectral weight in the pres- 
ence of a t' term in the band dispersion. Here W(D, T) 
(Eq. E3, dash-dotted line), W DDW (D, T) (Eq. gZJ (solid 
line) are shown together with W(T) (dashed line) which is 
for the non interacting case. The parameters are t' = 0.3 1, 
Tddw = 0.12 1, D(0) = AT DDW , and the doping 5 = 0.1. 
The chemical potential is evaluated self consistently at each 
temperature. Observe that near Tddw a small decrease of 
W DDW with respect to W(T) is seen, due to the chan ge of 
chemical potential near Tddw- Adapted from Ref. Il2(j 



to the Hamiltonian (|39[) . where Ak is the superconduct- 
ing order parameter. Ak = Ao [cos(fc x a) — cos(k y a)] with 
Ao the superconducting gap amplitude. In this case the 
OS in units of 7re 2 /h 2 is given by 



W UUW (D,A,T) 



„ MBZ 

-Y E U 

N V 



- - k tanh 



— — tanh — — 



2T 



(49) 
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FIG. 29: (Color online) Optical spectral weight in units of 
e 2 n/h 2 in the normal state (dash double-dotted line), in 
the DDW state Eq. (solid line) and in DDW+SC state 
Eq. I49H (dashed line). The values of parameters are for the 
doping S = 0.13 \ D (T = 0) = 0.92 1, A(T = 0) = 0.064 1, 
see Refs. Il20l and Il2ll for further details. The critical tem- 
perature is marked by the arrow. Observe that the decrease 
of W DDW {D, A,T) below T c is small. Inset: spectral weight 
plotted as a f unction of (T/t) 2 . Adapted from Ref. Il27l (see 
also Ref. EH. 



with E±k = J ^ k 

In Fig. we show results reproduced from Benfatto 
and SharapoAiiS-i for the effect of combined DDW and 
superconducting transition described by Eq. The 
optical weight is given in units of Z (nearest neighbor 
hopping) as is the temperature. The dash-double-dotted 
curve [W(T)] is for reference and gives results for the 
tight binding band without interactions while the solid 
lines gives W DDW (D,T) as before and the dashed curve 
is for W DDW (D, A,T) which includes the effect of a d- 
wave superconducting gap of amplitude A. We note that, 
as expected, this leads to an increase in KE with respect 
to the pure DDW case and, therefore, a drop in the OS. 

It is clear that, as yet, a complete theory of the OS 
in the DDW model does not exist. The calculations of 
Aristov and Zeyher 130 with vertex corrections properly 
accounted for give the opposite temperature dependence 
than found experimentally. This theory, however, does 
not include correlations beyond those directly responsible 
for the DDW transition. It is clear from what we have 
described here that correlations leading to lifetime effects 
need to be included as these are very closely related to the 
observed temperature dependence of the OS and cannot 
be ignored. 



VI. OPTICAL SPECTRAL WEIGHT 
DISTRIBUTION 

In this review we focused mainly on the OS for a sin- 
gle band integrated over all energies of relevance. As 
is seen from Eq. this quantity can be computed 
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FIG. 30: Selection of conductivity spectra for three BSCCO 
samples, underdoped, T c = 70 K (B70KUND, top frame), op- 
timally doped, T c = 80 K (B80KOPT, middle frame), and 
overdoped, T c = 63 K (B63KOVR, bottom frame). For the 
B70KUND sample, within the showed spectral range, the 
spectra at 50 K and 10 K are indistinguishable. Adapted from 
Ref.Hl 

from a knowledge of the single electron spectral den- 
sity A(k, oj) which determines the probability of occu- 
pation nk.o- of the state |k, a). This is a much simpler 
problem than computing the frequency dependent con- 
ductivity from a Kubo formula which involves the two- 
particle Green's functions. However, calculating the full 
frequency dependent conductivity cannot be avoided if 
one wishes to discuss the optical spectral weight distri- 
bution. The partial integration of o\ (T, lo) to a maxi- 
mum lo equal to lom has proved be very useful and has 
provided valuable information about normal and super- 
conducting state beyond what is obtained from the OS 
itself. We give here only two examples. In Fig. [201 we 
reproduce the conductivity data of Santander-Syro et 
al& in three BSCCO samples, underdoped T c = 70 K 
(B70KUND), optimally doped at T c = 80 K (B80KOPT), 
and overdoped at T c = 63 K (B63KOVR) at the various 
temperatures noted in the figure. From these data it 
is possible to calculate W(T, %) = J Q+ M duo <ti(T, lo) at 
various temperatures. Here the + indicates that the 
delta function representing the condensate has been left 
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FIG. 31: Rato AVF/Ws vs frequency showing the exhaustion 
of the Ferrell-Glover-Tinkham sum rule at conventional ener- 
gies for the OVR (diamonds) and OPT (triangles) samples. 
An unconventional (~ 16000 cm -1 or 2eV) energy scale is re- 
quired for the UND sample (circles). Note that the frequency 
scale changes at 800 and 8000 cm -1 . The changes in spectral 
weight are taken between 80 K - 10 K, 91 K - 10 K, and 100 K 
- 10 K for the OVR, OPT, and UND samples, respectively. 
Adapted from Ref. 15a. 



out. In Fig. 1311 we reproduce the experimental results 
of Santander-Syro et al£& for the normalized change in 
spectral weight AW/Ws vs lom in 10 3 cm -1 based on the 
data of Fig. |20| Open diamonds, triangles, and circles are 
for B63KOVR, B80KOPT, and B70KUND, respectively. 
Note the two breaks in the horizontal scale for lom at 800 
and 8000 cm -1 . The change in optical spectral weight 
is slightly different for the OVR, OPT, and UND sam- 
ples as indicated in the caption. Note the approach to 
its asymptotic value of one. For the overdoped case the 
scale is ~ 600 cm" 1 while for the UND sample it is much 
larger and of order 2 eV. These experiments clearly reveal 
a fundamental difference in behavior between overdoped 
and underdoped samples. This was also observed in the 
YBCO series. Data from Homes et alM* are reproduced 
in Fig. for electric field E parallel to the a-axis in 
YBCOe.gs (solid line) and YBCO 6 . 60 (dotted line). For 
the optimally doped sample the OS is rapidly saturated 
[u>m ~ 800 cm -1 ) while for the underdoped case a fre- 
quency of about 9000 cm -1 is needed. 

For a BCS s-wave superconductor the expectation is 
that the saturation of the OS should occur at a frequency 
lo a few times the gap A even if the system is dirty with 
scattering rates ~ 2A, Refs. Il32l and Il33l However, su- 
perconductors are better described by Eliashberg theory 
which properly accounts for coupling of the electrons to 
phonons. In this case the weight in the coherent quasipar- 
ticle part of the spectral function is Z = 1/(1 + A) where 
A is the mass enhancement factor. The rest of the spec- 
tral weight lies in an incoherent phonon induced band at 
higher energy, usually in the infrared. This part of the 
spectral function A(k, lo) contributes the so called Hol- 
stein band to the optical conductivity. Only the quasi- 
particle part is included in BCS theory, yet for A > 1, say, 
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FIG. 32: (Color online) The normalized weight of the conden- 
sate [N n {u)) — N s (u))]/ps,a vs lu for optimally doped YBCO6.95 
(solid line) and underdoped YBCO6.60 (dotted line) along the 
a-axis direction. The condensate for the optimally doped ma- 
terial has saturated by ~ 800 cm -1 , while in the underdoped 
material the condensate is roughly 80% formed by this fre- 
quency, but the other 20% is not recovered until much higher 
frequencies. The error bars on the curve for the underdoped 
material indicate the uncertainty associated with the Ferrell- 
Glover-Tinkham s um r ule. Inset: The low-frequency region. 
Adapted from Ref. Il3l 



more of the electron spectral weight is in the incoherent 
part than one finds in the coherent part. This part intro- 
duces a new energy scale into the problem, namely, an 
average phonon energy and it is no longer true to say that 
readjustment of optical spectral weight on entering the 
superconducting state can only occur on the scale of twice 
the gap. In fact, one should expect that when the elec- 
trons pair, an absorption process involving a phonon uje 
would be changed and shifted to an energy of 2 A + uje, 
thus shifting the Holstein band to higher energies and, 
hence, the incoherent part of the optical spectral weight 
in the superconducting state extends to higher energies. 
This implies that Jg+dv [ctn(T, v) — o~s(T, v)\ for large uj 
should saturate from above rather than from below as is 
the case in BCS. This is known from the s-wave phonon 
mediated casein and is illustrated in Fig. 1331 taken from 
Carbotte and Schachingeai^ for a <i-wave superconduc- 
tor using for I 2 x(oj) an MMP form. What is shown in 
the top frame are numerical results for W(T,u) defined 
as W(T, uj) = f a+ dv <j\{T,v). The real part of the con- 
ductivity <7i(T, uj) in arbitrary units on which these var- 
ious curves are based are shown in the bottom frame. 
These calculations are all done for infinite bands and 
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FIG. 33: Top frame: Optical spectral weight W(uj,T) = 
Jq + dv ai (v) for various cases as a function of uj. The dot- 
ted (dash-double-dotted) curve is for the normal state at 
T = 10 K (T = 95K), the solid curve for the superconduct- 
ing state at T = 10 K. The dashed (dash-dotted) curve is the 
difference curve between superconducting at T = 10 K and 
normal states at T = 10 K (T = 95 K). The approach of the 
difference curve to its saturated large u) value depends sig- 
nificantly on the temperature used for the subtracted normal 
state. The thin dash-double-dotted horizontal line is the value 
of the condensate contribution (penetration depth). The bot- 
tom frame shows the real part of the conductivity for the 
normal state at T — 293 K (dashed curve), T — 95 K (dash- 
dotted curve), T — 10 K (dotted curve), and for the super- 
conducting state at a T — 10 K (solid curve). All curves are 
for YFJCO6.95 with the impurity parameter given in Ref. 



corresponding Eliashberg equations for a <i-wave super- 
conductor so that the Ferrell-Glover-Tinkham sum rule 
holds, i.e.: the total optical spectral weight is conserved 
between normal and superconducting state^Z^ Param- 
eters were varied to get a good fit to data on YBCOg.gs- 
The reader is referred to the paper of Schachinger and 
Carbotte^ for details. Some impurity scattering in the 
unitary limit is included to get Fig. 123 W(T, uj) is shown 
for uj up to 250 meV. The solid curve is in the supercon- 
ducting state at T = 10 K and the dotted the normal 
state at the same temperature. Wn(T) (normal state) 
rises much more rapidly at small lu than does WsiT) 
(supercond. state) and goes to much larger values. The 
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difference W N (T = 10K,w) - W S (T = 10K,w) (dashed 
curve) is the amount of optical spectral weight that has 
been transferred to the condensate between (0 + , u>). This 
curve rapidly grows to a value slightly below the hori- 
zontal line representing the condensate contribution to 
the total sum rule. After this the remaining variation is 
small with a shallow minimum around 30 meV followed 
by a broad peak around 60 meV which falls above the 
thin dash-double-dotted horizontal line before gradually 
falling again towards its asymptotic value which must be 
equal to the condensate contribution. All these features 
can be understood from a consideration of the curves for 
a i (T, oj) given in the bottom frame. Comparing dotted 
and solid curves, we see that they cross at three places 
on the frequency axis at u>i ~ 8 meV, u>2 ~ 32 meV, and 
u>3 « 130 meV. These features are a result of the shift 
in incoherent background towards higher energies due to 
the opening of the superconducting gap. 

In an actual experiment it is usually not possible to ac- 
cess the normal state at T = 10 K (say) so that Wn(T = 
95K,w) just above T c = 92 K needs to be used (dash- 
double dotted curve). The difference Wn(T = 95 K, oS) — 
Ws(T = 10 K, oj) is shown as the dash-dotted curve in 
the top frame and is seen to merge with the dashed curve 
only at higher values of ui. Reference to the bottom frame 
shows that the real part of the normal state conductiv- 
ity at T = 95 K (dash-dotted curve) is much broader 
than at T — 10 K (dotted curve) and this accounts for 
the slower rise towards saturation of the dash-dotted as 
compared to the dashed curve in the top frame. We note, 
however, that W N (T = 95K,w) - W S (T = 10K,u) still 
approaches the penetration depth curve from above but 
much of the structure seen in the dashed curve is lost by 
using the data for the normal state at T = 95 K rather 
than at T = 10 K. Nevertheless, the energy scale over 
which the condensate is formed is set by the energy of the 
spin fluctuation spectra which extends up to 400 meV in 
our calculations and not by the gap value 2Aq. 



VII. SUMMARY 

The temperature dependence of the OS has recently 
been measured up to ~ 300 K in the normal state of 
several cuprates as well as in the superconducting state. 
Some controversy remains about details such as the ex- 
act temperature dependence followed in the normal state. 
However, it seems established that underdoped samples 
show anomalous behavior when they enter the supercon- 
ducting state. For a normal BCS superconductor the KE 
should increase while it is seen to decrease over its ex- 
trapolated normal state value. This, on its own, does not 
mean that the condensation is purely driven by the ki- 
netic energy but it does mean that we are faced with 
non BCS behavior. For optimally doped samples the 
change in KE is close to zero with a definite crossover 
to conventional behavior in overdoped samples. For the 
normal state many experiments, but not all, give a T 2 



dependence of the OS. Recent DMFT calculations for 
the Hubbard model provide evidence for a universal T 2 
behavior in the normal state. We argue, however, that 
such a law is not robust when electron-boson theories for 
the interaction are considered. In this case, if the boson 
energy is low, a linear in T dependence can result and 
by implication other dependences on T can arise for dif- 
ferent electron-boson spectral densities. The same holds 
for coupling to spin fluctuations in the NAFFL model. 
The issue of the temperature dependence of the under- 
lying normal state OS or KE impacts the analysis of the 
change in KE that results from the superconducting con- 
densation. If the normal state KE decreases faster at 
low temperatures than it does above T c (as we have seen 
in some models) this could be interpreted as an anoma- 
lous superconducting state, yet it would not be. On the 
other hand, in the NAFFL model as well as in electronic 
models more generally, anomalous behavior of the opti- 
cal sum can be understood as due to the so called col- 
lapse of the inelastic scattering rates in the superconduct- 
ing state if this reduces the effect of spin fluctuations at 
small lu. Such a hardening of the spin fluctuation spec- 
trum is thought to be a viable explanation for the peak in 
the microwave conductivity which is observed at a tem- 
perature considerably below T c . This effect, which can 
be modeled by a low energy cutoff in the electron-spin 
fluctuation spectral density results directly in a decrease 
in KE of the underlying normal state as additional co- 
herence sets in with the establishment of superconduc- 
tivity. This idea is related to the work of Norman and 
Pepin 87,88 who model directly the charge carrier self en- 
ergy rather than go through the electron-boson spectral 
function. They take the imaginary part of the self energy 
(£2) to be large in the normal state (incoherence) and in- 
crease coherence in the superconducting state through a 
low frequency cutoff below which £2 is zero. Another 
related model, which is even simpler, is to use a temper- 
ature dependent elastic scattering time which collapses 
in the superconducting state as aT 4 with a taken from 
microwave experiments. All these mechanisms can eas- 
ily provide savings in KE sufficiently large to explain the 
OS experiments. None, however, are fully quantitative at 
present. They also need to be extended to the overdoped 
side of the phase diagram. In this case one would expect 
the corresponding microwave peak to be greatly reduced 
and be more conventional. More sophisticated, but per- 
haps less transparent numerical calculations based on the 
Hubbard model and on the t-J model have also found 
that, as a function of doping, the KE can favor pairing 
in underdoped systems while for the overdoped case it 
behaves as for a conventional superconductor. Earlier 
calculations based on the Hubbard model, however, give 
kinetic energy pairing even on the overdoped side of the 
phase diagram. By contrast, calculations based on the 
negative U Hubbard model used to describe the BCS-BE 
crossover do capture the observed change in sign of the 
kinetic energy difference. 

In Ref. Il33l and, more recently, in Ref. the possi- 
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bility is raised that some of the effects studied in this 
review may arise from infinite bands when a finite cutoff 
is applied to the optical sum integral as it must in the 
analysis of experimental data. Certainly, for the high T c 
oxides there is no clear 'ending' of one band at a given 
energy followed by a gap before the next band sets in at 
higher energies. Instead, only a minimum is observed in 
the real part of the conductivity as a function of lu at 
about 1.2 eV. This energy has been taken as the cutoff 
on the transitions associated with the single band of in- 
terest. This cutoff is also roughly consistent with what is 
known about the width in energy of the bands in the ox- 
ides. Nevertheless, ambiguity remains and a theoretical 
study of overlapping bands and what this might imply for 
the single band sum rule would be valuable in clarifying 
this situation further. 

Much remains to be done to connect weak coupling ap- 
proaches with the more numerical strong coupling results. 
There is also a need to get more accurate data which 
could resolve the experimental debate that still goes on 
and to achieve quantitative agreement with experiments. 
What is clear, however, is that such experiments can give 
useful information on the correlation effects that are in- 
volved and they certainly have shown that superconduc- 
tivity on the underdoped side of the phase diagram of 
the cuprates is likely to be unconventional. It cannot 
be understood on the basis of a simple BCS model or an 
Eliashberg model with conventional phonons as these are 
not expected to change much as superconductivity sets in 
and they are distributed over a large energy range rather 
than peaked at low to. However, we stress again that as 
we have seen in one case, a small upturn in the OS at 
low temperature in the superconducting state can arise 



in spin fluctuation models without a readjustment of the 
electron-spin fluctuation spectrum, when the spin fluctu- 
ation energy in the MMP spin susceptibility is sufficiently 
small. While this upturn is less in the superconducting 
state than in its normal state at this same temperature 
(below T c ) it is above the value extrapolated to T = 
from the normal state above T c . A second important con- 
clusion of this review is that lifetime effects are important 
in determining the temperature dependence of the OS in 
both normal and superconducting state. These have to 
be accounted for in any definitive theory of this effect. 

There is a need to continue to explore other models 
such as phase fluctuations. In particular, the -D-density 
wave model as a possible competing order for the pseu- 
dogap phase has not yet included effects of correlations 
beyond those that lead to the DDW order, i.e: lifetime 
effects which we have argued to be central to this prob- 
lem. More work involving other models of interaction 
effects so as to understand what other temperature laws 
are possible for the OS vs T would also be useful. 
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